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This paper considers the motion of test bodies described in a simplified manne- - 
by means of multipole moments of the energy-momentum tensor density. Such a der 
scription was first employed by Mathisson in 1937. In the present paper the Mathisson 
method has been considerably simplified as a result of the use of the Dirac Ó function 
in the energy-momentum tensor. A new invariant definition of the centre of mass is 
also given. This definition ensures the uniqueness of the equations of motion. 


1. Introduction 


The general problem of the motion of bodies in general relativity theory is very 
difficult. Only approximate solutions to this problem in rather simple cases are known. 
This is the reason why one often considers simplified problems based on the investi- 
gation of the motion of the body of a.mass small enough in relation to the masses of the 
surrounding bodies so that the effect of the presence of this body on the gravitational field 
is negligible. We shall call such a body a test body. An example of test bodies are the 
planets in the field of the Sun. A further simplification of the problem is the description 
of test body by means of multipole moments of the energy-momentum tensor density, 
where only the most essential moments are taken into account. The simplest example 
of a multipole particle is the point test particle. The equations of motion of test 
particle were derived from the equations of the gravitational field by Infeld and Schild 
(1949). The general method of deriving the equations of motion of multipole particles 
in the linearized theory of the gravitational field was given by Mathisson (1931). 
A far less perfect (noncovariant) method was used by Papapetrou (1951) in general 
relativity theory. In the present paper we shall draw on Mathisson's method, simplifying 
it and extending it. A completely new approach to the problem of the motion of multi- 
pole particles will be the subject of a forthcoming paper. We shall consider here only 
multipole particles representing macroscopic bodies. Similar methods were employed 
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to obtain a classical model of elementary particles. One should mention here, above all, 
the papers of Frenkel (1926), Hénl and Papapetrou (1939, 1940), Mathisson (1940) 
and Weyssenhoff (1947). 


2. Equations of motion of test bodies. Definition of a multipole particle 


From the equations of the gravitational field 
G” = — 82 O” QU 


where G”’ is the density of the Einstein tensor and 0” the tensor density of the energy- 
-momentum of the matter giving rise to the field, there follow, by the Bianchi identities 


p, Ge 0 (2.2) 
four equations 


p, oO” —0 (2.3) 


which are the integrability conditions of the field equations. 

Equations (2.3), together with the equations of thermodynamics and of the theory 
of elasticity — not resulting from the field equations — determine the motion of 
matter. We shall call them the equations of motion. Thus in the general theory of rela- 
tivity the equations of motion need not and cannot be postulated separately, since 
they are contained in the field equations. This fact is of great importance; at the 
same time, however, it complicates the problem of motion, since the field equations 
have a solution only if the equations of motion are satisfied. On the other hand, the 
equations of motion cannot be solved without knowledge of the field. This difficulty 
vanishes in the case of a test body with a mass negligible as compared to the mass 
of the surrounding bodies. Let us suppose that the energy-momentum density of 
a body is proportional to a certain arbitrarily small parameter e. In the limit e = 0 this 
body will be precisely the test body, sincè, in the limit, the gravitational field will depend 
only on the remaining bodies. Dividing the energy-momentum density by & we obtain 
a finite density O"" representing the geometrical structure of the test body. We shall 
call this test density the energy-momentum density, remembering that it is defined to 
the accuracy of a constant coefficient. For every finite e, O"" must satisfy the equa- 
tions of motion; these equations will also be satisfied for € = 0 despite the fact tha 
they cease to be integrability conditions of the field equations. ; 

The description of test bodies by means of the continuous energy-momentum 
density OT is not the only possible description. 

Another description of the bodies equivalent as regards the equations of motion 
can be obtained by means of multipole moments of the energy-momentum density. 
Let us take an arbitrary time-like world line x” = £" (t) running inside or in the vicinity 
of the world tube occupied by the body and replace equations (2.3) by an equivalent 
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system of equations 


fn, 9" dq x= 0 
2 (67) 


(2.4) 


Í Ger (x^ — Er) p, OM da x = 0 


where Q (ET) denotes the three-dimensional neighbourhood of the point x” = & (t) 
(for a fixed t) containing the body. 
Taking to account the fact that 


V, Ox Se 2, 0% an Tm 0%, (2.5) 


expanding the Christoffel symbol on the hypersurfaces ¢ = constant into a Taylor 
series about the point x” = EI 


~ — 1 
IL E 9, I. (X — £) + 51 CA DÀ, ten — £5) (af a +... 


HET CM M C ENS, 


and carrying out the integration, we obtain, instead of the system of equations (2.4), 
the equations 

d 
ip. a 0^ HO, STS pi + Be D psd t= 0 


Te Taper 
- pok EE St ga — qm DU + 9 2, DL, try = 
G 


E [nim o9 Ln 5 Em, tri T0) o6 — nti mt L- 


SE T [ri rn HA CR 8, rn, [non $xÀ NEE 0 

j ee fx ue ce m (2.7) 

where £^ 7n" are the multipole moments of the density O"" with respect to the point 

eilt): 

open — f (xh 8)... an Er) 9" dey x. (2.8) 
a (er) 

We thus see that in the equations Hi motions the multipole moments are fully 


sufficient to describe the body. The description by means of a full, infinite set of 
multipole moments called the gravitational skeleton is not at all simpler than the 
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description by means of the continuous energy-momentum density. In certain physical 
problems the description by means of a continuous density is however too detailed. 


In such cases we can use a description by means of multipole moments, keeping 
only the most essential of them and neglecting the remaining ones. Thus, for example. 
a small planet moving in the field of the Sun is normally treated as a point test particle. 
which is equivalent, as we shall show, with the neglecting of the multipole moments 
higher than the unipole. If the rotation of the planet is taken into account the dipole 
moments should be retained and the departure from spherical symmetry can be 
described by means of quadrupole moments, etc. A body whose gravitational skeleton 
contains only a finite number of selected multipole moments will be called a multipole 
particle. Such a particle does not exist in reality; it represents only the most essential 
characteristics of a real body in some physical problems. (The concept of multipole 
particle also does not correspond to any specific approximation method in the sense 
of an expansion into a series with respect.to a small parametr). 

The equations of motion of a multipole particle can be obtained from Eq. (2.7) 
by retaining in them a finite number of expressions. The method outlined above was 
used by Papapetrou (1951). The calculations involved in this method can be somewhat 
facilitated by applying Eq. (2.4) to the singular energy-momentum density 


Tes — ef à — web 9, Ô + E PhD 20 Es (2.9) 


instead of the continuous density @%. Breaking off the series (2.9) would correspond 
to passing over to a multipole particle. The most serious shortcoming of the Papapetrou 
method is its non-covariance. 

One of the methods of avoiding this shortcoming is to describe the multipole 
particle by means of a clearly covariant singular energy-momentum density 


co 


T -f E GB — pr, (VP 0) +... + 


— CO 


c 


Virus tte Ög) [a (2.10) 
and at the same time replacing Eq. (2.4) by the invariant equation 


A Rg CU (2.11) 
D 


The multipole moments of the density (2.10), beginning with the order n + 1 
vanish in any arbitrary system; indeed this aul therefore represents a multipole 
particle. m i 

A description of multipole particles by means of the densis (2.10) corresponds 
to the method of Mathisson (1937), the shortcoming of this description being the lack 
of direct interpretation of the coefficients tı--Yx=# (corresponding to the coefficients 
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me? in Mathisson’s paper). In the case of unipole and pole-dipole particles, which 
we shall discuss later on, the coefficients t appearing in the final equations of motion 
have a simple interpretation. The more complicated multipole particles require, 
however, another description. In a forthcoming paper we will introduce a better 
description of the multipole particles and apply it to the pole-dipole-quadrupole 
particle. 

In the next section we shall give some auxiliary theorems concerning the singular 
tensor densities of the type (2.10). With the help of these theorems we then derive from 
Eq. (2.11) the covariant equations of motion of unipole and pole-dipole particles. 


3. Auxiliary theorems 


A. The singular tensor density 


n 


Ne m — > f Vu. EE (5) An (x? EE (3.1) 


k=0 


with arbitrary coefficients af: can be reduced to a canonical form 


È J ess tete (5) Beg G^ — POI d. Së 


where mir zm satisfy the conditions of symetry 
per Apr Hm m gU) be 2 e (3.3) 
and orthogonality to the vector Es 
gi Hi Km Gm. = 0, k= il. Boo Ue (3.4) 


B. The equations i à 
f Nim KR. A da x = 0. (3.5) 
D 


where K, ym is an arbitrary tensor field and D is an arbitrary four-dimensional 
region, lead to the vanishing of all coefficients of the canonical form of the density 
Nm. 

The importance of these theorems for the equations of motion of multipole 
particles becomes obvious if one notes that the expression V, T™ has precisely the 
form (3.1). Reducing this expression canonical form we obtain, by theorem B, the 
equations 9f motion. 

Theorems A and B were formulated in considerably more general form than is 
nesessary to obtain the equations of motion. These theorems can, however, also apply 
to other cases, e. g. in electrodynamics to the continuity equations 

j | 


pP =0 | ne a (8:6) 
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if the current j^ is represented as a singular vector density : 


M= 


je i V». WS) Ae (x — EI (s))] ds. (3.7) 


k=0 


We shall prove theorem A in the next section for those cases in which it will 
be necessary to obtain the equations of motion of unipole and pole-dipole particles. 
The proof will be based mainly on the formula: 


oo od 


| V n fae (s) & Are (x4 — &^ (s))] ds = Í — ôw ds (3.8) 


LÉ — 00 


whose correctness may readily be verified in the coordinate system which is geodesic 
along the world line x* = & (s). 

The proof of theorem B is immediate. Inserting (3.2) into (3.5) and integrating 
we obtain : 


$2 


n pen 
21 f 1) V... OR phe em ds = 0 (3.9) 
Gab Y m 


where s, and s, correspond to the points of the world line lying on the boundary of D. 


We pass over to the rest system (£, = 0, Ey = 1) in which, by (3.4) we have 
n Ss — 
Diem eds). (3.10) 
Ko H k m 


Making use of the arbitrariness of K»,...«,, we assume that all space derivatives 
of the order lower than n vanish. We now obtain 


Ge k K yi itm ds = 0. à (3.11) 


A 21.40% 


$3 " 
J m 
5 


. From the arbitrariness of the nth space derivatives it now follows that 
prety rey — () (3.12) 
in the rest system, or 


plan Ma Km () (3.13) 


in any arbitrary system. ; 
The vanishing of the remaining coefficients v is established in a similar manner. 


4. Equations of motion 
` a) Unipole particle 
A test body of vanishing dimensions is called a test particle. The most proper 


description of such a particle by means of multipole moments is obtained by taking 
as the world line of the particle the line x” = £'(s). The only non-vanishing multipole 
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moments in this case will be the unipole moments; one should therefore ascribe to 
the test particle the energy-momentum density 


T” = J v" (s) O | v" = c" (4.1) 


— C0 


which means that the test particle is simply a unipole particle. 
Let us write 7^ in the form 


ET = p) gu E | m” È L m’ Er L m" (4.2) 


where 
m = 1“ D e E at E Er E) £, > 
mi” — el (8 — E E) (6% — E E). (4.3) 
m” and m“ are obviously orthogonal to vector Et, 
In: order to reduce the expression 


prc i V, l(m£ & + m’ £ + m" E + m”) ôg] ds (4.4) 


— 09 


to canonical form it is sufficient to make use of Eq. (3.8). We obtain at once 


co 


p,T" = ne (m + m^) dea) + y» Wan + m’ Er) dca] | ds. (4.5) 


— 00 


Eq. (2.11) therefore gives, by theorem B, the following equations: 
m" + m” gH =) 


Sa. 
ES (m&# + ml = 0 (4.6) 
‘that is 
mi e mn =O 
dm Aën : | 
np E —— = 0. 4.7 
ds E de : eu 
Collecting our results, we obtain for the test particle: 
T. = f m gr E da de 


m = constant 
ay (4.8) 


The world line of the test particle is therefore a geodesic line. 
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BM eee 
The constant m can be interpreted as the rest mass, since in the rest system 
m= JT dx. (4.9) 

QE) 


The mass of the test particle vanishes in accordance with the definition of this 
particle for uncharged particle it has no meaning. If we examine a charged particle 
with a test charge e, then the ratio e/m does not vanish and the magnitude of this 
ratio determines the motion. 


b) Pole-dipole particle 


In accordance with the definition the pole-dipole particle is described by the 
density 


oo 


T= [f [s G^ — £ ()) — DT Se) (^ — EI ()1. 


T” a= T, Tr = rt (4.10) 


as in the case of a unipole particle we represent t” in the form (4.2). A similar decom- 
position of t“” gives in the general case 


v" — n* EM fr Ln" EL am ër ne Ex q” (4.11) 
where 
nt = 77 EE 
am — 700 (8t — EXE) (OE — E £) &,, 
nw" — e (dz — EE) (0, Er E) (5 EE 
sch. (4.12) 
Employing theorem A we can set q” = 0, the coefficients 7"" being changed correspond- 
ingly. In order to simplify the interpretation later on, it is better to take 
ger = ` WEEN (4.13) 


which is always possible by (3.8). 
Introducing the notation 


S" = mn” + wb E nt EH (4.14) 


we now have 


1 1 
qum — neu A. = S" Er + = Sx Eu, (4.15) 
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In reducing p, T” to canonical form we make use of the results for a unipole 
particle. The remaining calculations are given in the Appendix. 
The following result is obtained: 


1 Agen . jae RS , 
parm = [IS Sn, = a EE CM og 


= : H de v Ey $ 
b 5 22 Re, | Ney a ZI + më” e e e 2) A a 


2 zb 


L pyy Er + n2») self ds. (4.16) 
We thus obtain the equations 
nore 4 s» EH = 0. (4.17) 


Multiplying (4.17) by E, we obtain 


ne = 0 (4.18) 
from which we have 
no =O (4.19) 
or 
ni = () (4.20) 
since 
p^ — n (4.21) 
The other equations 
f 1 Agen aco d 
v SCENE 8 — E Ey 4.22 
me mr Er (Ei) (4.22) 


can be split, as previously, into a part parallel 


120577 
TAE 4.23 
m^ = 23875 = , ( ) 
and orthogonal 
oi SL Nis 
m" — E a (85 — È” En) (04 — & EA = 0 (4.24) 


to the vector €. The tensor S% is antisymmetric by (4.14) and (4.18). 
Separating Eq. (4.24) into symmetric and antisymmetric parts we now obtain 


E (4.25) 
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and 
I ine) Ob uin) 
ds 
= a + & ës et x (4.26) 
From the results obtained thus far we have further 
E n? | = J | s S9 R" ua & = 0. (4.27) 


We use Eqs. (4.20), (4.23), and (4.25) to reduce the number of unknown quantities 
characterizing the particle. Taking into account these equations we write the energy- 
-momentum density in the form 


co 


SES Ke 1 a KK d às SHA 
tem EE EE Tr Su s) de 


— 3 VA Sh & + s de dëi m (4.28) 


The quantities m and S appearing here have a simple interpretation. A simple 
calculation gives 


Ji — mn — tw — Py To) dg s — 3: (4.29) 
and 
fest a, =" (4.30) 
in any arbitrary system, and 
TRE (4.31) 


in a rest system. The quantities 5", S" and m should therefore be interpreted as the 
internal angular momentum, statical moment, ahd rest mass of the particle, respecti- 
vely. The remaining equation, that is Eqs. (4.26) and (4.27) are identical with the 
equations given by Papapetrou (1951). 

The physical state of the particle is described by the values of the mass, internal 
angular momentum, statical moment, and the derivative of the statical moment 
representing the position of the centre of mass of the particle with respect to the world 
line. The equations of motion should have a unique solution for any given initial 
state, and therefore should contain the first derivatives of the mass and angular mo- 
mentum and the second derivative of the statical moment, which indeed occurs, as 
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readily verified, in the rest system of the world line in which the equations of motion 
take the form 


ke 
E 0 
à? Sro ~ 1} d 
ds? ar mnl, + DJ S% Roxa =0 
d AS 
TF + ml = 0. (4.32) 


The world line with respect to which the multipole moments are counted have 
been regarded thus far as completely arbitrary. We can draw conclusions about the 
motion of the centre of mass only from knowledge of the mass and the statical moment. 
A considerably clearer picture of the motion could be obtained by relating the line 
x” = E" (s) to the centre of mass. The world line of the centre of mass of a particle in 
a given system will be the line with respect to which the statical moment vanishes. 


Eqs. (4.26) and (4.27) should therefore be supplemented by the condition 
Sl (4.33) 


and treated as equations for the mass, angular momentum, and coordinates of the 
centre of mass Pr. Condition (4.33) was used by Corinaldesi and Papapetrou (1951) 
who investigated the motion of pole-dipole particles in a Schwarzschild field. The 
coordinate system to which this condition referred was the rest system of the Schwarz- 
schild singularity. The equations of motion with condition (4.33) have a unique solution 
for any arbitrarily given initial state determined ty the value of mass, angular momen- 
tum, and coordinates and velocity of the centre of mass. Condition (4.33), which 
is simultaneously a definition of the centre of mass, can be written in covariant form 


S" K, =0 (4.34) 


where K, is a field of time vectors. This condition depends on the coordinate system 
or on the field K,. In particular, one may require that condition (4.33) be satisfied 


in the rest system of the world line, which corresponds to taking &, for the vector KS 
S" E, = 0. (4.35) 


Condition (4.35) could serve as an invariant definition of the centre of mass 
if it would be single-valued. It turns out, however, that for each particle there exists an 
entire class of world lines which satisfy (4.35) and the equations of motion. The non- 
-uniqueness of condition (4.35) is conected with the fact that the vector &,, in contrast 
to the vector K,, cannot be given in advance, since it is the quantity we are seeking. 
Mathisson (1939) who introduced the equations of motion of pole-dipole particles 
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pov eU EE AR I E EA RE 


under the assumption n* = 0 in formula (4.11) obtained the equations 


E E 
Sx t ‚Sur » — () 
ds > ds a ds d 
ër Sur u 1 “A RE. ey; 
m gs SU ds 5 Se ee) 
Si D =0, m= constant. (4.36) 


There are indeed as many of these equations as there are unknown quantities. 
The non-uniqueness is contained in the order of these equations, which is too high. 
Besides Mathisson, condition (4.35) was used by Weyssenhoff (1947) and Pirani 
(1956, 1957). The non-uniqueness of this condition can be illustrated by a simple 
example which will serve for the construction of an invariant and unique definition 
of the centre of mass. 

In the case of flat world the general equations (4.26) and (4.27) permit as one 
of the solutions the solution 


een eo 
fe ee 


GO Se B, — 0 i = Constant (4.37) 


We shall show that any solution in the flat world can be obtained from this solution 
by a suitable change of the world line. Passing to a Lorentz system, in which the 
particle is at rest, we write the solution (4.37) in the form 


= == 0,0 EIOR a N A o (4.38) 


The energy-momentum density representing the particle has in this system the 
form 


T = m ô (x) 
T"? = : Srs 2, dra (at) 


T° =0. ` (4.39) 


The same particle described by means of multipole moments with respect to 
an arbitrary world line x“ = n” (s) has the energy-momentum density 


T = mg (x — n) + ma 8, Aa (X — qf) 
T — Z S79, Ne — 1) | 
T" = 0, (4.40) 


Denoting the rest mass, internal angular momentum, and statical moment of the 
particle with respect to the line x“ = ^ (s) by m, S" and S” respectively, 
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we obtain 


dt 


De pe 


Sr = Sr, SO = my (4.41) 


where ds’? = g,, dn“ dn”. We can verify that the equations of motion are now satisfied 


automatically by (4.38). For the new values in a flat world and Lorentz frame! Eq. 
(4.26) gives 


dSrs Së dS§s4 s. dS?r 
REUS DEER a UE 


= mn n°? ns mm wenig 0 (4.42) 


and 


dSr0 dS% dëi 
^r Li ^0 , 
ds’ = N N ds? + NN ds’ 


en mn" Ets, mn” ne 95 + m? n° UM = 


= m7" (Sw v^ n" — 1) = 0. (4.43) 
Similarly, Eq. (4.27) gives 
d Gs Au dsr D d / H H H . 
i (mnt + E n) mtem = 0 (4.44) 
and 
duse s d5” f d mu Passt 
Se (ore mu) ege mr = 0 (4.45) 


The calculations above show that the world line taken to describe the particle 
is indeed entirely arbitrary. We see, further, that the conditions 


SK, = S® K,— my K,=0 
S” E, = mm K,—0 d (4.46) 


uniquely determine the world line for given S” and K,, the second of these conditions 
being a consequence of the first. On the other hand, the conditions 


m 


Sn, = — S" n° —mnn’=0 


S” yn, = — mn nt =0 (4.47) 


allow of an entire class of world lines. In order to determine this class we divide the 


dnt 


1 , 
n " denotes " 
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dt od z 
first of these conditions by gz dy and rewrite in vector motation: 
e mos GE 
HO AR (4.48) 


where S (S38, S91, $1?) is the vector of the internal angular momentum, and r (oi, 92,19) 
is the radius vector. 

It follows from (4.48) that the radius vector 7 is perpendicular to the velocity 
vector, i. e. it has a constant length. 3 

In addition, the radius vector is perpendicular to the fixed vector S, and the motion 
takes place along a circle in the plane perpendicular to S, where between the radius 
of the circle r and the velocity v on the circumference we have the relation 


mro SV (4.49) 
which means that the angular velocity of the motion is 


m 
UE c (4.50) 

The class of solutions obtained here corresponds to the solutions of Weyssenhoff 
(1947) derived directly from the equations of motion. 

Let us return to the problem of the invariant definition of the centre of mass. 
Condition (4.34) with the field not connected with the motion of the particle is not an 
invariant, and the requirement that K, — n, allows of an entire class of world lines 
of which only one is the world line of the true centre of mass. 

Among the lines satisfying the condition 


Sn = 0 (4.51) 


one is singled out; it is the geodesic line x“ = £" (s). It is natural to regard precisely 
this line as the world line of the true centre of mass. 
This line satisfies the condition 


ST? = 0 | (4.52) 
which can be written in an arbitrary system in the form 


S^ p, = 0 (4.53) 
where iq 
Ó Sv 
ds 


is the momentum vector of the centre of mass. Condition (4.53) is both invariant and 
single-valued. 


pr =m = mn’ + —— mh (4.54) 
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Returning to the usual notation, we shall show that Eqs. (4.26) and (4.27) in 
flat space with the condition 


So (4.55) 
where 
Lp ee 
p = mi + ba (4.56) 


have a unique solution. In flat space Eq. (4.27) can be written in the form 


äer _ 
A =o (4.57) 


Multiplying Eq. (4.26) by p, we obtain immediately 


ISA 
= | 
CS &,=0 (4.58) 
from which we get 
p-—m&. (4.59) 
The equations of motion now take the simple form 
6S” Aë : 
Up Un atus 0, SwE,—0, m = constant. (4.60) 


It is characteristic of these equations that there is a lack of any coupling between 
the progressive motion of the particle and the angular momentum. In flat space there 
may occur only an apparent coupling between the progressive motion and the angular 
momentum. The coupling of this type appearing in Mathisson's equations (4.36) is 
partially apparent, and partially real, the former being connected with the nonunique- 
ness of condition (4.35) and the latter, with the curvature of space. In this case the 
equations of motion in general do not allow of solutions in the form of a geodesic 
line. Owing to the lack of another definition we keep the definition (4.55) also in 
a curved space and we shall investigate the form taken by the equations of motion 
of a centre of mass defined in this way. Instead of (4.57) we now have 


1 f 
P Lis Ra B= 0. (4.61) 


Instead of (4.59) we obtain in a similar way 
2 AES À ers 
mi = TE SA (6 — e) 


€ 


+ = Sa suo Ryan (8% — & &) &. (4.62) 


Use of (4.62) in the equations of motion complicates rather than ‘simplifies them. 
Let us now suppose that the particle moves in a Schwarzschild field and that the order 
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of the linear dimensions R of the particle is much less than the order of the distance 
from the singularity L. The right-hand side of (4.62) is in this case proportioned to 
the square of a small parametr e = R/L. Other terms of this order would appear in 
the equations of motion in which the quadrupole moment is taken into account. Since 
we have neglected terms connected with the quadrupole moment we should, to be 
consistent, also neglect the right-hand side ef equation (4.62). As the equations of mo- 
tion we now obtain the equations 
KE 

"de 


fest 
+ 54 Rh P0 


m = constant 
Sr sect 


óSw" 
ds 


TONO SA . + SA 
P Gn EA see Sa em 0. (4.63) 


Further simplification of the last of these equations would lead to contradic- 
tion with the previous ones. 

Equations (4.63) are the result of employing the multipole precedure along with 
an approximation method. This would be justified not only in the case of a Schwarz- 
schild field, since the description of the particle by means of the unipole nad dipole 
moments is correct only for the cases in which the right-hand side of Eq. (4.62) is 
actually small. 


Appendix. 


Let us reduce to canonical form the expression 
x En Eu 1 v £x Eu 1 u Ex En xv E au £v 7 
V vx n EE ee E D EL ns Er 1 p dc ds (A.1) 


where the coefficients n are orthogonal, in all indices, to the vector Er. Using Eq. (3.8) 


we obtain 
à Ô T 2 1 è 
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In the second brackets we change the order of differentiation 
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1 Ó she = 
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pen Prom esc 4 m9 da d (A3) 
We again apply Eq. (3.8) to the last brackets 
^ Sc . : Get 
i G br DA we Pp am Pr A ntm) Rs te + 
i OSE ; 
Sc i V» tds. (a ar Vx (ES Eu + n*"") 2 ds. (A.4) 


We represent the derivative of the antisymmetric part of the last brackets by 
means of the curvature tensor 
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o) H Eu [(00 + n) EH) dal) ds (A.5) 


Making use of the identity 
[vx] a 0 (A.6) 


in the first brackets and applying Eq. (3.8) to the second we arrive at the final result 
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VERHALTEN DES OFFENEN LUFT-SPITZENZAHLERS BEI HO- 
HEREN TEMPERATUREN 


Von T. Lewowski unD B. SUJAK 


lustitut für Experimentalphysik der Universität Wroclaw 
(Eingegangen am 20 Dezember 1958) 


Ein heizbarer offener Versuchs-Luftspitzenzähler wurde konstruiert. Das Verhalten 
des Zählers wurde bei Temperaturen aus dem Bereich 20°C — 500°C studiert. Es wurde 
festgestellt, dass auf einfachem thermischen Wege die Arbeitsspannung des offenen 
Luftspitzenzählers um ca 30% erniedrigt werden kann. 

Der Luftspitzenzáhler wurde auch im Bereich der Dauerentladung betrieben und 
sein Verhalten bei den Temperaturen von 20?C bis 500?C untersucht. Es wurde fest- 
gestellt, dass der koronierende Spitzenzühler (Funkenzähler) der auf a-Teilchen anspricht, 
bei Temperaturen über 100°C vollkommen unabhingig vom Wasserdampfinhalt der 
Luft ist. Der Zahler zeigt dabei ein flaches und langes Plateau auf seiner Arbeitscharakte- 
ristik und eine recht kleine Arbeitsspannung, die bei ca 450°C in der Nahe von 2300 V 
liegen kann. 


Einführung 


Die offenen Luft-Zähler, die als unstabil arbeitende Detektoren bekannt sind, 
werden bei den Kernphysikalischen Uritersuchungen nicht gern benutzt, besonders 
nach der Entwicklung der Szintillationstechnik. Jedoch in den letzten 10 Jahren sind 
diese offene Luft-Zahler wiederum der Gegenstand mehrfacher Untersuchungen ge- 
worden. Diese haben unter anderen auch zur Konstruktion einfacher &-Detektoren 
geführt (siehe z. B. Simpson 1948; Frank, Lynck und Baumgardner 1955; Burek, 
Lewowski und Sujak 1957). Auch die Probleme der Registrierung von Exoelektronen 
haben auf die Notwendigkeit der Erweiterung der Kenntnisse über den offenen 
Zähler, besonders über den offenen Luft-Spitzenzähler hingewiesen (Kramer 1949, 
1950; Bohun 1955; Sujak 1957, 1959). 

Aus einigen Arbeiten (siehe z. B. Simpson 1948, Frank und Mitarbeitern 1955) 
ist bekannt, dass eine kleine Erhóhung der Temperatur des Zählers über die Raum- 
temperatur sich sehr günstigt auf den Hintergrund des offenen Luft-Zählers auswirkt. 
Die Messungen der Thermostimulation-Coelektronenemission haben aber auf die 
Änderungen der Arbeitsspannung des sich während der Messung erwärmenden 
Spitzenzählers hingewiesen (siehe z. B. Sujak 1958). i 

Die Temperatur erweist sich also als ein wichtiger Faktor der Arbeitsbedingungen 
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auch des offenen Luft-Zählers. Dieses veranlagte uns dazu das Verhalten des offenen 
Luft-Spitzenzihlers gegenüber einer Erhöhung der Temperatur näher zu untersuchen 
und der vorliegende Bericht enthält die Ergebnisse dieser Untersuchung. Dabei 
wurden besonders zwei Ziele verfolgt: 

1) Die Möglichkeit die hohe Arbeitsspannung des offenen Luft-Spitzenzählers 
auf dem thermischen Wege zu erniedrigen 

2) Das Verhalten des offenen Luft-Spitzenzählers zu studieren, der in der soge- 
nannten Dauerenladung betrieben wird (Funkenzähler) und als einfacher «-Teilchen- 
Detektor arbeitet. 


Versuchsbedingungen 


Der benutzte Versuchs-Luftspitzenzihler wurde in einen kleinen elektrischen 
Widerstandsofen eingebaut. Der Zylinder des Zählers von 18 mm Durchmesser wurde 
aus Silberblech angefertigt. Als Spitze wurde ein Platindraht von 0,5 mm Durchmesser 
benutzt, der am Ende zu einer Halbkugel abgerundet oder zu einem Kegel verschärft 
wurde. Um die Erwärmung des Isolators (Polystyren), in dem die Spitzenhalterung 
eingebetet wurde, zu verhindern wurde sowohl die Spitzenhalterung als auch die 


Ia 8 7 


Abb. 1 Schematische Darstellung der- Versuchsanordnung: 1-Millivoltmeter, 2-Autotransformator, 

3-Hochspannungsspeisegerat, 4-Verstürker, 5-Untersetzer, 6-Impulsmittelmesser, 7-Polistyrenpropfen, 

8-Zählerspitzenhalterung aus Neusilber, 10-Zahlerspitze aus Platindraht, 11-Zählergehäuse aus Silber- 
blech, 12-Thermospitze des Thermoelementes. 


Zylinderverlängerung wegen der schlechten Wärmetransporteigenschaft aus Neusilber 
angefertigt. 

Die Temperatur des Versuchs-Zählers konnte im Bereich von 20°C bis 500°C mit 
Hilfe Autotransformators geändert werden. Sie wurde mit Hilfe eines Cu-Konstantan- 
Thermoelementes gemessen. Als Temperatur des Zählers wurde die Temperatur 
angenommen, die mit Hilfe des gennannten Thermoelementes gemessen wurde. Die 
Lage der Thermospitze des Thermoelementes sowie die ganze Versuchsanordnung ist 
schematisch auf Abb. 1 wiedergegeben. 

Der elektronische Teil der Versuchsanordnung bestand aus dem elektronisch 
stabilisierten Hochspannungsspeisegerät, dem Proportionalverstärker mit regulier- 
barem Verstärkungsfaktor, dem Untersetzer und Impulsmittelmesser (Abb. 1). 
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Arbeitsspannung des Zählers als Funktion der Temperatur 


Wie bekannt, haben die Untersuchungen der Durchschlagsspannung in Gasen 
unter anderen zu der Aufstellung des Gesetzes von Paschen geführt (siehe z. B. Loeb 
1950 Seite 367). Laut diesem Gesetz ist die Durchschlagsspannung V, im Gas eine 
Funktion des Gasdruckes p und der Elektrodenentfernung d 


Vg = f(p:d) (1) 
Diese Funktion stellt beinahe eine Proportionalität für Luft dar, wenn Vz die 


Werte von ca 1000—8000 V annimmt (Loeb 1950 Seite 368). Für diesen Fall kónnen 


wir also schreiben 
EC pud (2) 
Bei einem elektrischen Durchschlag hängt der Einfluss des p: d von der Zahl 
der Moleküle die einem Elektron auf jeden cm des Weges zukommt. Dieses hängt 


aber bei einem Konstantem Druck von der Dichte D des Gases ab. Wir kónnen also 
nach Loeb (1950 Seite 514) in die Proportionalität (2) anstatt p: d, D: d einsetzen 


Vg—D:d bei p = const. .... (3) 


Wenn p sich nicht ändert, kann die Dichte D des Gases in erster Náherung als propor- 
tional zu 1/T angesehen werden. Wir kónnen also erwarten, dass die Durchschlagsspan- 
nung reziprok der Temperatur geündert werden kann: 


Vg — E bei p — const (4) 


Wenn wir das Problem des elektrischen Durchschlages näher betrachten, so ist 
es durchaus zu erwarten, dass auch die Proportionalitat (4) der Einsatzspannung für 
selbstständige Entladung (Dauerkorona) in erster Näherung zukommt. Vor dem 
Bereich der selbststindigen Entladung (der Dauerkorona) liegt aber der Bereich der 
unselbststindigen Entladung, den eigentlich aber der Spannungsbereich des Geige- 
rzühlers darstellt. Wir kónnen also eine recht grosse Herabsetzung der Arbeitsspannung 
V, eines offenen Luftspitzenzählers (nach Geiger) auf rein thermischem Wege zu 
erwarten und zwar reziprok mit der Temperatur. 

Um den Einfluss der Erwürmung auf die Arbeitsspannung des Luft-Spitzen- 
zühlers zu verfolgen, wurden bei den jeweiligen Temperaturen die Arbeitscharakteris- 
tiken des Záhlers aufgenommen. Ein Radium-Priparat diente als Quelle jonisierender 
Teilchen bei der Aufnahme von Arbeitscharakteristiken des Versuchszühlers. In der 
Strahlung waren also sowohl &-wie ß-Teilchen wie auch y-Quanten vorhanden. Die 
Geometrie des Zählers sowie die Lage des Radiumpräparates wurden während einer 
Reihe von Messungen festgehalten. i 

Soweit die Genauigkeit der Spannungsmessungen erlaubte, wurde bis über 450°C 
festgestellt, dass der Charakter der Arbeitscharakteristik des Zählers nur wenig von 
der Temperatur abhängt. Im wesentlichen zeigen die Charakteristiken nur eine 
Verschiebung in der Richtung der kleineren Arbeitsspannungen. Bei zirka 480°C 
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Abb. 2 Einfluss der Temperatur auf die Arbeitscharakteristik des Versuchszählers. 
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beginnt jedoch eine Thermoemission schon die Arbeitsweise des Zählers zu stören. Die 
Störung beruht auf einem beträchtlichen Anwachsen des Hintergrundes des Zählers. 
Beispiele von 2 aufgenommenen Arbeitscharakteristiken für 40°C und 330°C sind auf 
Abb. 2 wierdegegeben. Der Wert der Anliegenden Spannung, der auf der Arbeits- 
charakteristik über der Einsatzspannung V, um etwa 25 V liegt, wurde als Arbeitsspan- 
nung V4 angenommen. Auf Abb. 3 ist die so definierte Arbeitsspannung V} als 
Funktion der Temperatur T für zwei Spitzen von verschiedenen Krümmungsradien 
wiedergegeben. Die Krümmungsradien der Spitzen wurden unter dem Mikroskop 
geschätzt, 

Es stellte sich heraus, dass laut der Erwartung, auf rein thermischen Wege sich 
durchaus die Arbeitsspannung eines offenen Luft-Spitzenzählers um etwa 1000 V 
erniedrigen lässt. Das bedeutet im Mittel eine 309 — tige Arbeitsspannungsenkung. 
Bei entsprechend dünnen Spitzen bedeutet das, dass stabilisierte Hochspannungsgeräte 
bis 2000 V zum Speisen eines „heissen“ Luft-Spitzenzählers ausreichen würden. 


Einfluss der Temperatur auf die Arbeitscharakteristik des &-Detektors 


Es ist bekannt, dass ein offener Spitzenzühler der sich in der sogenannten Dauerent- 
ladung befindet, die einzelnen &-Teilchen registrieren kann, indem Streamers gebildet 
werden. Dabei werden grosse Stórrungen notiert, die dem in der Luft sich befindenden 
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Abb. 4 Temperaturabhingigkeit des „Hintergrundes“ des koronierenden Versuchszühlers. Das Profil 
der Pt-Zählerspitze ist schematisch wiedergegeben. Ohne a-Teilchen-Quelle. Relative Feuchtigkeit der 
Luft etwa 40%. 
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Wasserdampf zuzuschreiben sind (Anderson und Hertz 1955, Anderson 1958, Burek, 
Lewowski und Sujak 1957, Sujak 1958 a, b). Der koronierende Spitzenzähler spricht 
auch mit Leuchtfäden auf Wasserdampf, sowie CO, in der Luft an*. Dieses Ansprech- 
vermögen ist selektiv gegenüber der anliegenden Spannung (Abb. 4), und ist als eine 
Störung zu betrachten. Diese Störung ist desto unangenehmer, dass sie sich nicht 
additiv gegenüber den Streamerimpulsen der a-Teilchen verhält. Es kann also nicht von 
einem wahren Hintergrund gesprochen werden, der dem Wasserdampf zuzuschreiben 
wäre. Die Werte des scheinbaren Hintergrundes lassen sich nämlich nicht von den 
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Abb. 5 Temperaturabhüngigkeit der Arbeitscharakteristik des koronierenden Versuchszühlers. Das 
Profil der eingesetzten Pt-Zählerspitze ist eingezeichnet. Bei allen Kurven ist die Lage der a-Teilchen- 
Quelle nicht geändert worden. Relative Feuchtigkeit der Luft etwa 40%. 


Messwerten einfach abziehen, wie das bei den Geiger-Müller Zählrohren üblich ist. 
Die Impulsenzahl, die dem Wasserdampf zuzuschreiben ist, ändert sich während der 
Messung von a-Teilchen sehr und lässt sich nicht kontrolieren (siehe z. B. Abb. 6). 

In einem vorgegangenem Bericht (Burek, Lewowski und Sujak 1957) haben wir 
bereits beschrieben, wie diese Störung einfach umgegangen werden sein kann. Der Zä- 
hler muss bei einer Spannung betrieben werden, die mehr als zirka 400—500 V über 
der Einsatzspannung für die Dauerentladung (Dauerkorona) liegt. 

Das Ansprechvermógen der koronierenden Spitze wurde inzwischen eingehender 


* Wie letzt durchgeführten Versuche zeigen spricht der koronierende Spitzenzühler auf me- 
hrere Flüssigkeitsdämfe an. Ein ausführlicher Bericht darüber ist in Vorbereitung (Anmerkung bei 
der Korrektur) 
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von Hertz und Anderson (1955, 1958) studiert. Aus den Berichten dieser Verfasser 
ist z. B. bekannt, das eine zum Kegel verschärfte Zählerspitze (Platin) besonders auf 
den Wasserdampfinhalt in der Luft empfindlich ist. Aus diesem Grunde eben haben 
wir in den Versuchsspitzenzähler auch eine scharfe Spitze eingesetzt. Wie unsere 
Untersuchungen gezeigt haben, wird bei der Erwärmung des offenen Zählers nicht 
nur die Einsatzspannung für die Dauerentladung (Dauerkorona) merklich herabgesetzt, 
aber es verschwindet auch das selektive Ansprechvermögen auf die Feuchtigkeit der 


Mt ' 
[Imp|sec] eeo 400°C 
— 25% 
2000 
3380 V 
7000 
0 


TR unde nt CR ES TR en). 


Abb. 6 Temperaturabhangigkeit der Zählrate N/t des Versuchszählers bei Anderungen der Entfernung 
d der a-Teilchen-Quelle vom Zähler. Die sonstigen geometrischen Verhältnisse wurden konstant gehalten. 
Die jewielige Temperaturen und Arbeitsspannungen sind bei den Kurven vermerkt. 


Atmosphärischen Luft (Abb. 4 und Abb. 5). Die unangenehme Störung bei der 
Benutzung eines koronierenden Luft-Spitzenzählers als a-Detektor kann also auf 
einfachem Wege durch Erwärmung des Zahlers beseitigt werden. Dabei kann die 
Arbeitsspannung des Spitzenzählers auf zirka 2300—2600 V herabgesetzt werden, 
wie das die Abb. 5 illustriert. > 


Abschliessende Bemerkungen 


Das durchgeführte Studium über das Verhalten eines offenen Luft-Spitzenzählers 
bei den Temperaturen 20°—500°C hat eindeutig gezeigt, dass auf einfachem thermi- 
schen Wege sich die Arbeitsspannung um zirka 30% herabsetzen lässt. Die theoretischen 
Erwägungen wurden experimentell geprüft und es konnte annüherend eine Herab- 
setzung der Arbeitsspannung des Zählers laut der Proportionalität V4 ~1/T fest- 
gestellt werden. Dies scheint auf den ersten Blick ein bischen unerwartet zu sein, da bis 
jetzt bei den Studien über Temperatureinfluss auf die Arbeitsweise der G—M Zähler 
immer über eine Erhóhung der Arbeitsspannung mit der Temperatur berichtet 
wurde. Man muss jedoch unterscheiden zwischen einem abgeschlossenen Zähler, wo 
bei einer Erhóhung der Temperatur der Druck des Zählgases wächst und damit auch 
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die Arbeitsspannung, und einem offenen Zähler wo der Druck des Zählgases konstant 
bleibt und die Dichte des Gases mit der Temperatur kleiner wird, was zu einer 
Herabsetzung der Arbeitsspannung führt. 

Durch einfaches Erwärmen des Zählers bis über etwa 100°C wird der Einfluss 
des Wasserdampfes in der atmosphärischen Luft neutralisiert. Dadurch erreicht man 
ein langes und flaches Plateau auf der Arbeitscharakteristik des &-Teilchen-Funken- 
zühlers mit einer recht kleinen Arbeitsspannung. Der Einsatz eines warmen Zählers 
ist besonders aussichtsvoll überalle dort wo sich eine höhere Temperatur der Arbeits- 
bedingungen nicht vermeiden lässt (z. B. Kernreaktoren). 

Es ist durchaus denkbar, dass sich ähnliche warme Luftfunkenzähler von grös- 
seren Oberflächen für «-Teilchen oder Neutronen (siehe z. B. Janik und Szkatula 
1956) herstellen lassen. Dabei könnten vielleicht die Zähldrähte selbst zugleich als 
Widerstandsheizkörper dienen. Es ist aber klar, dass dabei einige technische 
Schwierigkeiten zu überwinden sind. 
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ZUR VORGETÄUSCHTEN THERMOSTIMULIERTEN COELEKTRO- 
NENEMISSION VON HYDRATEN 


Von B. SuJAK*) 


Institut f. Physik d. Polnischen A. d. W., Wroclaw; Institut f. Experimentalphysik d. Universität Wroclaw; 


und 


A. BOHUN 
Institut f. technische Physik d. Tschechoslovakischen A. d. W., Prag. 
(Eingegangen 23 Dezember 1958) 


Es werden „Wasser-Glow-Kurven“ für einige Hydrate, wie Na,SO, : 10 H,0, 
BaCl,2 H,O, Na,S,0,:5 H;O und CuCl,2 H,O mit Hilfe eines im Koronagebiet arbei- 
tenden Luftspitzenzühlers mit scharfer Spitze gemessen. Gleichzeitig wird der Verlauf 
der Temperatur des Heizofens photographisch registriert. Es konnte gezeigt werden, 
dass beim benutzten Heizofen von kleiner Wärmekapazität der Temperaturverlauf unter 
den Maxima der ,,Wasser-Glow-Kurve“ nicht linear ist. Dieses weist auf die endother- 
mische Zersetzungsreaktionen der untersuchten Hydraten hin. 

Es wird gezeigt, dass auch ein Geiger-Zühler mit dem Kügelchen auf den durch 
ein Hydrat abgegebenen Wasserdampf indirekt mit Impulsen anspricht. 


Einführung 


Sujak (1958 a, b, c) hat gezeigt, dass die Abgabe von Kristallisationswasser 
wührend der thermischen Zersetzung der Hydrate, die auf thermostimulierte Coelektro- 
nenemission untersucht werden, eine Exoelektronenemission leicht vortäuschen kann. 
Der offene Luftspitzenzähler, der sehr oft bei dem Studium der thermostimulierten 
Coexoelektronenemission benutzt wird, spricht nümlich auf das durch die Probe 
freigegebene Kristallisationswasser mit Maxima der Zahlrate an 

Wenn im Zähler mit scharfer Spitze bereits Koronaentladung herrscht (so genannte 
Dauerentladung), so spricht er mit „Streamers“ nur noch auf alfa-Teilchen und Was- 
serdampf an, und es kónnen keine einzelne Exoelektronen mehr registriert werden. 
Mit Hilfe eines auf solche Weise arbeitenden Zühlers mit scharfer Spitze kónnen doch 
Eh >: 032 2 E eer a E 

*) Der Verfasser ist dem Tschechoslovakischen, sowie dem Polnischen Ministerium für Hoch- 
schulwesen für die finanzielle Unterstützung seines Aufenthaltes in Prag, wo diese Arbeit gemacht wurde, 
sehr verbunden. 
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die „Wasser-Glow-Kurven“ analog zu den Exoelektronen-Glow-Kurven gemessen 
werden (Sujak 1958 a, b, c). 

In dem vorliegenden Bericht wollen wir kurz über weitere Versuchsergebnisse 
bei den unbestrahlten oder mit Röntgenstrahlen bestrahlten Proben berichten, die 
mit Hilfe einer anderen Apparatur im Institut f. technische Physik d. Tsch. A. d. W 
erziehlt worden sind. Bei diesen Messreihen konnte gleichzeitig mit dem Verlauf aer 
vorgetäuschten Emission auch der Verlauf der Temperatur registriert werden, was 
bei den früheren Messungen nicht der Fall wai. 


Versuchsbedingungen 


Die Messungen wurden mit Hilfe der Apparatur durchgeführt, die von Bohun 
(1955) beschrieben wurde im Zusammenhang mit den Untersuchungen von Halogeni- 
den der Exoelektronenemission. Die Zählrate des Zählers, sowie der Verlauf der 
Temperatur wurden gleichzeitig photographisch registriert. Die elektronischen Kreise, 
die zum Messen der Lumineszenz bei den Exoelektronen-Untersuchungen dienten, 
wurden ausgeschaltet. Die Temperatur des Heizofens wurde automatisch mit einer 
Geschwindigkeit von ca 2 grad/sec erhóht. Der Durchmesser des Luftspitzenzählers 
betrug 2 cm. Das offene Ende des Zählers wurde mit einem Metallgitter abgeschirmt. 
Die Maschendichte des Gitters betrug etwa 150 Maschen/cm?. 

In dem Geiger-Zähler, der gewöhnlich mit einem Kügelchen von zirka 0,2 mm 
Durchmesser beim atmosphürischen Luftdruck (Spannung etwa 2100 V) arbeitet, 
wurde dieses durch eine gewöhnliche Grammophonnadel ausgewechsellt. Die Arbeits- 
spannung des Zählers wurde um etwa 300 V über der Einsatzspannung für Dauer- 
entladung festgelegt, und der Zähler sprach nicht mehr auf die f-und y-Strahlung an. 
Die Impulse der Koronalavinen wurden mit Hilfe eines Diskriminators abgeschnitten 
und nur die Impulse der Streamers wurden durchgelassen und registriert. 

Die kleinen unbestrahlten, (Orientierungsweise auch mit Röntgenstrahlen 
bestrahlten) Kristallstücke der Hydrate wurden unter soeinem Zähler mit Hilfe des 
Heizofens gleichmässig erwärmt und die Temperatur des Ofens sowie die Zählrate des 
Zählers Ny,/t dauernd registriert. 

Alle Messreihen wurden im Dunkeln durchgeführt. 


Versuchsergebnisse 


Die Messungen wurden an folgenden Hydraten durchgeführt: Na,SO,-10H,0, 
BaCl,2H,0, Na,S,035H,O und CuCl,:2H;0. 

Auf Abb. l ist der zeitliche Verlauf der Temperatur beim Ofen ohne Probe 
wiedergegeben. Der Verlauf ist beinahe linear. Als aber auf den Heizofen ein kleines 
Kristallstiick Na,SO,-10H,O gelegt worden ist, so wurde mit dem Luftspitzenzähler 
mit scharfer Spitze auf der Wasser-Glow-Kurve ein ausgepregtes Maximum in der 
Nahe von 350 K einregistriert (Abb. 2). Für Wasser-Glow-Kurve pflegen wir ebenfals 
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die Bezeichnung Ny,/t = f(T) zu benutzen (Sujak 1958 b). Bei dem Temperaturverlauf 
ist nun im Gebiete von Maximum sehr starke Unregelmessigkeit zu verzeichnen, die 
für eine endothermische Reaktion spricht. Obwohl die Wärmeenergie, die in jeder 
Zeiteinheit dem Ofen zugeführt wird, die gleiche ist wie für den Fall in Abb. 1, zeigt 
die Temperatur manchmal sogar Absinkentendenzen. Dies bedeutet, dass die Wärme- 
kapazität des Ofens verhältnissmässig klein ist und deswegen die nötige Wärme zum 
Ausdampfen des Kristallisationswassers sofort eine Anderung der Temperatur zur 
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Abb. 1 Zeitlicher Verlauf der Temperatur des Heizofens (T) und des Hintergrundes des Zählers (A). 


Folge hat. Wenn nun dieselbe Probe zum zweiten mal erwärmt wird, tritt auf der 
Wasser-Glow-Kurve schon kein Maximum mehr hervor und der Verlauf der Tempe- 
raturkurve ist wiederum linear. 

Ein Beispiel der Messreihen für Na,S:035H,0 ist auf Abb. 3 wiedergegeben. 
Wenn die Entfernung zwischen dem Zähler und der Probe (d. h. Gitter-Probe) kleiner 
gewühlt wird, so erscheinen die entsprechenden Maxima breiter und oft bei kleineren 
Temperaturen. Weitere Beispiele der Wasser-Glow-Kurven zeigen die Abb. 4 und 
Abb. 5. Abb. 4 gibt eine typische Wasser-Glow Kurve für BaCl,:2H,0 und Abb. 5 
für CuCl,-2H,O an. ER, 

Um festzustellen, ob ein Geiger-Zähler, dessen Anode mit einem Platinkugelchen 
vom zirka 0,2 mm Durchmesser versehen ist, tatsächlich auf das abgegebene Kristal- 
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400° 
Abb. 2 Wasser-Glow-Kurve: Ny/t = f(T) für Na,SO410 H,O. Zähler mit scharfer Spitze. Entfernung 
Gitter-Probe etwa 3 mm. 
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300? 850? 400° 450° 500°550°K Zeit 
Abb. 3 Wasser- KC are Ny/t = f(T) fiir Na,S,03°5H,0. Zähler mit scharfer Spitze. Entfernung 
à Gitter-Probe etwa 3 mm. 
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Abb. 4 Wasser-Glow-Kurve: Ny,/t = f(T) für BaCl,-2H,O. Zähler mit scharfer Spitze. Entfernung 


Gitter-Probe etwa 2 mm. | 
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Abb. 5 Wasser-Glow-Kurve: Nyy/t = f(T) für CuCl,2H;0. Zähler mit scharfer Spitze. Entfernung 
Gitter-Probe etwa 3 mm. 
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lisationswasser anspricht, wurde folgender Versuch angestellt: der Zähler wurde mit 
gewóhnlicher Spannung von etwa 2100 V betrieben. Er sprach also auf die B-und y- 
Strahlung an und konnte die einzelnen Exoelektronen registrieren. Nun wurde ein 
kleines Kristallstück von Na,S,03'5H,O auf den Heizofen gelegt und eine normale 
thermostimulierte Coelektronenemissionskurve aufgenommen. Diese hat ein scharfes 
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300° 350? 400°450390° 550°K Zeit 
Abb. 6 Die vorgetäuschte Exoelektronenemission Glow-Kurve: Nit = F(T) für Na,S,03°5 H,O. Zähler 
mit einem Kügelchen. Entfernung Gitter-Probe etwa 1,5 mm. 


Maximum in der Umgebung von 450°K (Abb. 6). Die Temperaturkurve wies dabei 
eine Unregelmässigkeit auf, die sehr eindeutig den sich abspielenden Prozess des 


Hydratzerfalles wiedergibt und für eine vorgetäuschte Exoelektronenemission-Glow- 
Kurve spricht. 


Abschliessende Bermerkungen 


Wie aus den durchgeführten Messreihen hervorgeht, wurden die früheren Fest- 
stellungen von Sujak (1958 a, b, c) auch mit Hilfe einer anderen Apparatur bestätigt. 
Es hat sich erwiesen, dass die Registrierung des Temperaturverlaufes sehr behilflich 
sein kann beim Unterscheiden einer vorgetäuschten thermostimulierten Coelek- 
tronenemission von der wahren Exoelektronenemission. Die besonderen Messungen 
von Wasser-Glow-Kurven können bei den Zweifelfällen zu entscheidenden Ergebnissen 


Zur Coeletronenmission von Hydraten ' 425 


fiihren, besonders dann, wenn keine photographische Registration des Temperatur- 
verlaufes vorgenommen wird. Es scheint, dass der lineare Verlauf der Temperaturkurve 
die Täuschung ausschalten kann. Das ist aber nur dann möglich, wenn ein Heizofen 
von sehr kleiner Wärmekapazität gegenüber der Wärmekapazität der untersuchten 
Probe benutzt wird, wie das z. B. im Falle der Apparatur von Bohun ist. 

Die Autoren danken Frau K. Listonovà für die Durchführung der Messreihen. 
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L'INFLUENCE DU RAYONNEMENT INFRAROUGE SUR LE 
TRANSFERT D’ENERGIE DANS ZnS-Mn 
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Institut de Physique de l’Académie Polonaise des Sciences 
(Recu le 7 janvier 1959) 


On a étudié linfluence du rayonnement infrarouge sur la photo-luminescence 
de ZnS-Mn en poudre à la température de —190°C. On a constaté l'affaiblissement de 
la bande bleue de la luminescence et le renforcement de la bande orangée. 

Les résultats de l'étude de la dépendance de ces effets de: la concentration de l'acti- 
vateur, l'épaisseur de la couche de phosphore, la longueur d'onde et l'intensité du rayon- 
nement excitant ainsi que de l'irradiation infrarouge montrent que les rayons infrarouges 
influencent le mécanisme du transfert d'énergie par les électrons et les trous. 


A basses températures on observe dans le spectre de photoluminescence dn ZnS-Mn 
excité par un rayonnement de longueur d'onde À = 3650 A une bande bleue et une 
bande orangée. Des centres liés aux perturbations du réseau ZnS — selon Kröger 
(1955) ce sont des lacunes de cations de charge positive — sont responsables de la 
bande bleue, des centres liés au manganèse étant responsables de la bande orangée. 
L'irradiation par les rayons infraronges du phosphore pendant l'excitation change 
d'une facon considérable la répartition d'intensités dans ce spectre (Jaszezyn 1957); en 
effet il en résulte un affaiblissement fort de la luminesnescence bleue et un renforcement 
de l'orangée. Un tel effet du rayonnement infrarouge montre qu'il exerce une influ- 
ence substantielle sur le transfert d'énergie entre divers centres de luminescence 
dans ZnS-Mn. 

Le ion bivalent du manganése Mn est connu comme activateur de la luminescence 
de beaucoup de cristaux de composition et symétrie différents. Sous l'influence de 
l'ultraviolet proche et de la lumière visible seuls sont possibles des transitions inter- 
dites ayant lieu dans la couche électronique 3d: de l'état 6S aux états 4G, *D, 4P et 4F. 
L'absorption correspondante à ces transitions est faible et la luminescence du manganèse 
a lieu avant tout aux frais de l'énergie transmise par d'autres centres qui absorbent 
le rayonnement excitateur plus considérablement. Ce sont des centres liés ou au 
réseau du cristal fondamental, ou aux impuretés introduites. Ainsi on a constaté 
(Botden et Kröger 1948, Botden 1952, Dexter 1953) que dans les phosphates, les 
carbonates et les silicates les ions du cérium, du plomb et d'autres transmettent par 
résonance l'énergie d'excitation aux ions du manganése. 

| (427) 
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Shogeo Shionoya et Kaiichi Amano (1958) supposent que l'effet sensibilisateur 
du Pb sur la luminescence du manganése dans ZnS-Mn, Pb est lié lui aussi à la réso- 
nance. 

En cas, où les centres pouvant fournir l'énergie au manganèse sont lonisés par 
l'irradiation excitatrice, un autre mécanisme du transfert à part de la résonance devient 
possible, notamment par l'intermédiaire de libres porteurs. De nombreux faits 
montrent que précisément les centres ,bleus" mentionnés ci-dessus, présents aussi 
dans ZnS nori activé, sont ionisés par le rayonnement de Wood: la présence de la 
photo- -conductibilité, le caractére de la cinétique de la luminescence et l'augmentation 
du rendement correspondante à l'augmentation de l'intensité excitatrice et à la baisse 


de la température en témoignent. 


Par effet de l'absorption d'un quantum d'ultraviolet un électron est arraché d’un. 


centre bleu, où apparait alors un trou localisé. La recombinaison de ce trou avec 
'électron libre est accompagnée d'un rayonnement bleu. Néanmoins, avant cette 
recombinaison le trou peut être libéré, par exemple sous l'action de l'énergie thermique; 
celà est d'autant plus probable que la température est plus élevée et l'excitation plus 
faible. Le trou libre, bien sûr, peut se lo caliser sur un centre bleu non excité, (Vinokurov 
et Fock, 1956) ou aussi recombiner sans rayonnement avec un électron localisé dans 
un piége. C'est pourquoi, la libération des trous affaiblit la luminescence bleue. Les 
quanta d'infrarouge ainsi que l'énergie thermique libérent des trous des centres bleus; 
d’où l'effet d'extinction exercé par l'infrarouge sur la luminescence de ZnS à basses 
température (fig. 1). 

Quand des niveaux liés 4 d’autres activateurs par exemple des atomes ou des 
ions ,interstitiels^ du manganèse apparaissent dans la bande interdite du ZnS, 
les trous libérés peuvent se localiser sur eux, transmettant de telle facon l'excitation 
des centres bleus aux centres orangés. Néanmoins MnS et ZnS sont isomorphes, 
et le manganèse apparait dans ZnS essentiellement sous l'aspect d'ions bivalents Mnt 
remplaçant les ions Zn** dans les noeuds du réseau. 

L'ionisation des centres Mn** demande plus d'énergie que celle des Zn** 
(à cause de la participation des électrons d à la liaison), le niveau fondamental de 
ces centresest situé sous la bande interdite les trous libres ne peuvent donc pas se 
localiser sur lui directement et le processus de stimulation par eux du rayonnement 
de ces centres est plus compliqué. 

Króger (1940) fait remarquer que la présence de ces centres déforme le potentiel 
périodique du réseau; ces perturbations deviennent des lieus de la recombinaison 
de l'électron et du trou, et l'excés d'énergie de l'électron est transmis par réabsor- 
ption (en vertu des résultats des dernières années on préfère parler de résonance) 
aux ions Mn**, provoquant St les transitions mentionnées aux niveaux de Stark. 

Quand les perturbations du potentiel du réseau causées par Mn sont suffisamment 
importantes elles peuvent jouer le róle des piéges de phosphorescence. Il semble que 
la recombinaison d'un trou libre et d'un électron situé dans un piége P pourrait 
exciter le manganèse directement par la voie suivante 
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Mutt + e — Mn* 
Mont + et — (Mntt)* 
(Mnt?)* Mnt + hy 


Un processus du méme genre a été suggéré en résultat l'étude sur d'autres 
phosphores (Antonov Romanovsky et Trapesnikova 1956). 

Le retour des ions Mn** à l’état fondamental est accompagné d'une émission 
orangée. 

On pourrait supposer que l'infrarouge stimule en partie le transfert par réso- 
nance ou par radiation de l'énergie des centres bleus aux centres orangés (par exemple 
grace à son influence connue sur la largeur de la bande de la radiation émise 
(Jaszezyn 1958b). Néanmoins l'évaluation de l'importance de cette influence dans 
le cadre du mécanisme décrit est impossible, car nous ne connaissons pas le spectre 
de l'énergie transmise. D'autre part les résultats expérimentaux ont l'air d'affirmer 
que l'influence de l'infrarouge sur le spectre de la luminescence est liée avant tout 
à laliberation des trous des centres bleus excités. e 

Sur les fig. 2 et suivantes représentant les spectres à — 190°C on lit sur l'abscisse 
les fréquences (en 101 sec-1), sur l'ordonnée des nombres N ' proportionnels au nom- 
bre de quanta émis corespondant de à la fréquence donnée. Il en résulte que les 
surfaces comprises dans les parties correspondantes du spectre entre les courbes de 
luminescence sans et avec irradiation infrarouge sont proportionnelles au nombre 
de quanta éteints dans la bande bleue AN, et gagnés dans la bande orangée. AN, 

Puisque, comme on a mentionné plus haut, les trous libérés des centres bleus 
et „perdus“ pour la radiaton bleue peuvent ou « recombiner sans radiation, ou OI 
recombiner avec excitation des centres orangés, (l'accoroissement du AN,), donc 
le coefficient u = (AN,/AN,) va changer dans diverses conditions de l'expériment 
selon la part relative des cas œ et f. 

(La valeur de u dépend non seulement de l'efficacité de transmission de l'excita- 
tion, mais aussi du rapport des probabilités de transition radiative dans un centre 
orangé et bleu). : 

Les fig. 2 et 3 représentent les spectres de luminescence du ZnS-Mn/10-#) et 
de ZnS-Mn (10-8) sans et avec irradiation infrarouge. On voit que l'augmentation 
du la concentration a une influence sur le changement du spectre ‘de la luminescence 
provoqué par l'infrarouge, notamment il a pour conséquences: Le 

I) l'affaiblissement relatif de la bande bleue W = AN,JN, (où AN, = LN? EE E 
—N,|) plus fort 

II) l'efficacité du transfert u plus grande SC 


Mmao?) > UNa’) H 
Par suite du transfert de l'énergie des centres bleus aux centres de l'activateur, la 
conséquence I de l'augmentation de la concentration de l'activateur dans le phospho- 
re est qualitativement compréhensible indépendamment. du mécanisme de ce transfert. 
Par contre le résultat IT ne peut étre expliqué que sur la base du mécanisme dés trous. 
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L'augmentation de la concentration du Mn cause l'accroissement de cette partie du 
nombre total des trous libres qui recombinent à proximité de ces centres en les ex- 
citant d'autre côté la part de recombinaisons sans rayonnement décroit. 
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Les diagrammes représentés sur la fig. 4 montrent l'influence de la longuer 
d'onde du rayonnement infrarouge sur l'effet discuté. On a étudié le changement du 
spectre de la luminescence sous l'influence d'une irradiation monochromatique infra- 
rouge de diverses longueurs d'onde, mais d'énergie approximativement identique. 

La comparaison des diagrammes 4 fait remarquer que le coefficient u est supérieur 
pour les quanta de l'infrarouge à l'énergie supérieure qui, comme le montre la table 1, 
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Table 1 
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1% — L’intensité stationnaire de la luminescence excitée par l'ultraviolet (A, = 3650 À). 
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. . H : RA cc 
provoquent un éclat plus fort, donc amoindrissent l'occupation des pièges („communs“) 
(Levchine, Rigikov 1958). Ces piéges ne sont pas lies localement aux contre: 
Mn. Un trou libre recombine sans rayonnement avec électron situé dans un tel piége. 
Le décroissement de l'occupation de ces pieges rend plus rares les recombinaisons sans 
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rayonnement, par le méme augmente le nombre de quanta produits dans la bande 
orangée ainsi que dans la bleue et par la suite augmente aussi la valeur du coefficient u. 

Une analyse plus précise de l'effet de l'irradiation infrarouge sur le spectre de 
radiation de ZnS-Mn exige de tenir compte de la répartition énergétique des piéges de 
tenir électroniques ainsi que de l'influence éventuelle de l'irradiation sur la transition 
des électrons des niveaux d'excitation d'un centre bleu à la bande de conductibilité. 
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La fig. 5 représente le spectre de luminescence d'une couche du ZnS-Mn (1073) 
a) épaisse (~ 1 mm) b) mince (— 10 my). Les deux échantillons ont été exécutés d'un 
matériel identique a) en comprimant une quantité correspondante de poudre dans 
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un châssis et b) en déposant sur le verre une émulsion à l'eau et en séchant ensuite 
dans l'air à la température de ~ 90°C. 

Comme il est facile d’apercevoir: 

1) Sans irradiation infrarouge le rapport des intensités de la luminescence orangée 
et bleue est inférieur dans la couche mince que dans l'épaisse 


(NN. > UNN, b 


p 


2) Le rapport des changements provoqués dans des parties diverses du spectre 
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par l'infrarouge n'est pas le méme non plus dans la couche „a“ et dans la couche „b“: 
(ANSAN pa > NAN) 


1) et 2) montrent pareillement que la probabilité du transfert de l'énergie de la 
bande bleue à l'orangée a diminué. Peut-étre la facon différente de préparation de la 
couche ,,b“ est favorable à l'apparition dans les graines de poudre d'états superficiels 
capturant les trous. Ce probléme exige de nouvelles recherches. Néanmoins il a paru 
intéressant de présenter encore un exemple de la variabilité de u, caractéristique pour 
le mécanisme des trous. La diminution de l'intensité du transfert par résonance ou le 
décroissement de la réabsorption dans une couche mince conduirait à la diminution 
du nombre de quanta „obtenus“ dans la bande orangée ainsi que „perdus“ dans la 
bande bleue. 

Remarquons que qualitativement le caractére de la dépendance des effets étudiés 
de l'intensité du rayonnement excitateur et de l'intensité de l'irradiation infrarouge 
est d'accord avec les prévisions faites sur la base du mécanisme des trous. 

Pour le phosphore excité plus fortement l'affaiblissement du rayonnement bleu 
par l'infrarouge devient plus petit car en résultat de l'augmentation de la concentration 
des électrons libres il devient plus probable qu'un trou localisé dans un centre bleu 
recombine, avant qu'un quantum infrarouge ne le rejéte vers la bande fondamentale. 

Le fait que les effets produits par le rayonnement infrarouge (AN, ANo) 
grandissent plus lentement que l'intensité de l'infrarouge peut étre facilement 
compris qualitativement méme sans exécuter les calculs correspondants. Il suffit en 
effet de remarquer que le nombre de trous rejetés par le rayonnement infrarouge 
des centres „bleus“ dépend non seulement de l'intensité de l'infrarouge, mais aussi 
du degré d'excitation du phosphore (c'est-à-dire du nombre total de trous localisés 
sur les centres bleus) qui tombe, bien sûr, avec l'irradiation infrarouge. 

Quand la luminescence est excitée par J'ultraviolet plus lointaih qui provoque 
des transitions bande-bande et qui produit donc directement des trous libres, on 
constate que le röle des trous „liberes“ par l'irradiation infrarouge diminue, done 
l'influence de ce rayonnement sur le spectre tombe, comme le montre le tableau 2. 


à 


Table 2 


Zn S — Mn (10-*) —190°C 


Aexcit (A) 


1% = L'intensité stationnaire de la luminescence (la bande orangée sans irradiation) 


peak = L'intensité statiohnaire de la luminescence (la banda orangée avec irradiation.) 
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Par contre l'intensité de la bande orangée excitée par la lumiere visible de longueur 
d'onde À = 5460À ne change pas sous l'influence de l'irradiation infrarouge. Le 
rayonnement de longueur d'onde 4 — 5460 À correspond à l'absorption caractéristique 
du manganése. 

Rappelons pour terminer que l'influence de l'irradiation infrarouge ne se limite 
pas au changement des intensités et à l'élargissement des bandes de luminescence. 
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Fig. 7a Le spectre d'émission du ZnS-Mn (10-4) non irradié par l'infrarouge 


Comme le montre la fig. 7, l'analyse du spectre de luminescence de ZnS -Mn(10-4), après 
avoir admis la symétrie des bandes, nous permet de constater outre la bleue et l'orangée 
la présence de 2 autres bandes: une rouge et une verte. En rapport avec la discussion 
qui s'est engagée sur l'existence de ces bandes dans ZnS-Mn (Klick 1955) remarquons 
que sous l'influence du rayonnement infrarouge l'intensité de la bande verte tombe 
comme celle de la bleue, tandis que celle de la rouge monte. Il se peut que le rayonne- 
ment „rouge“ est émis par des centres de manganèse situés dans les positions inter- 


stitielles. , 
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A STATISTICAL MOLECULAR THEORY OF ELECTRIC, MAGNE- 
TIC AND OPTICAL SATURATION PHENOMENA IN ISOTROPIC 
DIELECTRIC AND DIAMAGNETIC MEDIA 


By S. KIELICH AND A. PIEKARA 
Institute of Physics, Polish Academy of Science; A. Mickiewicz University, Poznan 
(Received February 16, 1959) 


It is the aim of the present investigation to establish a unified statistical-molecular 
theory of the nine electric, magnetic and optical saturation phenomena in isotropic 
dielectric and diamagnetic media (gases, condensed gases, liquids). For the case 
of a condensed medium composed of polar molecules of arbitrary symmetry, anisotro- 
pically polarizable and non-linearly deformable in an external field, genera] expressions 
yielding the nine molar constants have been derived, namely: 
group I — the electric saturation: S% in an electric, Sy in a magnetic and S% in 

i an optical field; i 
group II — the magnetic saturation: Sy in an electric, sm in a magnetic and SY 
in an optical field; 

group III — the optical saturation: 55, in an electric, S?" in a magnetic and Sy in 
an optical field. 
A discussion of these molar constants is given for particular cases of spherical and 
axial symmetry of the molecules. For axial symmetry, the general formulas reduce to 
those given previously. Moreover, expressions have been derived which relate 
the above-mentioned molar constants to the variations of the ‘electric permittivity 
As‘, Ae”, Ae°, of the magnetic permeability Au‘, Au”, Ay’, and of the optical 
refractive index An“, An", An? of the medium respectively, as resulting from the 
action thereon of a strong polarizing electric, magnetic or optical field. Only three 
of the nine possible effects under consideration have been detected to-date, namely, 
the electrooptical Kerr effect, the magneto-optical Cotton-Mouton effect and the 
electric saturation in an electric field, i. e. electro-electric saturation. The authors 
have derived equations for computing each of the six as yet unknown quantites Je”, 
Ae’, ... An? from the known experimental Kerr and Cotton-Mouton constants. By 
these formulas, the quantities Ae”, Ae°, ... An? have been numerically computed 
for nitrobenzene. 

(439) 
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1. Introduction 


The electric, magnetic and optical properties of an isotropic medium are deter- 
mined by the three scalar coefficients e, u and n characteristic of the substance. 
The electric permittivity e and the magnetic permeability y of the medium are measured 
applying an electric or magnetic feld, respectively, which vary slowly, and whose 
frequency is much smaller than that of the Debye dispersion. The optical refractive 
index n is measured by an optical field, which is here understood to be the electric 
field of an electro-magnetic wave whose frequency exceeds that of the Debye dispersion. 
Such fields will be termed the electric E’, magnetic H^ and optical £ measuring field, 
respectively. 

The present paper deals exclusively with dielectric, diamagnetic and isotropic 
media. Placed within one of these fields, the medium undergoes polarization of a kind 
corresponding to the measuring field and. determined, respectively, by the vectors 
of electric polarization P*, magnetic polarisation P" and optical polarisation P". 
The relation between the electromagnetic material coefficients e, u and n, on the one 
hand and the vectors P*, P" and P? on the other, is given by the following general 
equations of electrodynamics: — 

e m 
e— 1= 4r se, DEE EE (1.1) 

Usually, measuring fields E’, H’ and & of small intensity are applied experimen- 
tally; as a result, each of the vectors P*, P" and P^ in eqs. (1.1) is a linear function 
of the respective measuring field giving rise to the polarization, and the quantities e, 
u and n are functions which do not depend on the field, but depend solely on the 
medium and its thermodynamically determined state. 

An entirely different picture is obtained if an additional external field of great 
intensity is applied, as e. g. a constant electric field E. In this case, the polarization 
arising from the external field tends to saturation; the polarization vectors P*, P" and 
P" are now generally non-linear functions of the external field, whereas the electro- 
magnetic material coefficients e, u and n now become functions (in at least the second 
degree) of the additionnally applied external electric field E. This latter field will 
be termed the polarizing electric field. 

The variations of the coefficients e, u and n as arising from the effect of a strong 
electric polarizing field E on the medium will be defined as follows: 


de — à ApS =o —p, Antes nn} (1.2) 


‚wherein e, u° and n° denote the electric permittivity, magnetic permeability and ` 
optical refractive index of the medium, respectively, as measured in the presence of 
the polarizing field E, whereas e, u and n denote the respective quantities as measured 
without a polarizing field (at E — 0). | 
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A constant magnetic field H, or an optical field 2, may serve, too, as polarizing 
field; hence, the variations of e, u and n may be defined in analogy to those of eqs. 
(1.2) as brought about by a magnetic polarizing field H: 


Ze =e" — AN = nu An" =n" —n, (1.3) 
or by a polarizing optical field £: 
Neve cct Ayo py AN nn. (1.4) 


Thus, a set of nine theoretically possible variations of the electromagnetic material 
coefficients €, u and n resulting from the effect of an external polarizing field E, H 
or & are obtained, as shown in the following table: 


polarizing field 


electric | magnetic | optical 


measuring field 


electric E’ 
magnetic H’ A ut 
optical €? À ne 


The variations defined in eqs. (1.2) — (1.4) will be termed, with respect to the 
coefficient measured, effects of electric Lie, As”, Jet, magnetic (Au‘, Au”, Au’) 
and optical (Anf, An”, An?) saturation. Strictly speaking, the effects Ae, Ae", 

. etc. found experimentally did not represent the true saturation effect because of 
two by-effects appearing in the measurements. These are the reversible quadratic 
electro- or magneto-caloric and electro-or magnetostrictive effects. By computing 
both by-effects (see Piekara, Chelkowski and Kielich, 1957) and subtracting from the 
effect measured, e. g. Jet, the true saturation effect is obtained, e. g., 


Aeg. = Ae! — A — À Eejectrostriek: (1.5) 


sat electrocal. 


In the case of electric saturation, the by-effects are found to be small, of opposite 
sign, and nearly cancelling. 


à 


I. Effects of electric saturation 


1. Electric saturation in an electric field, i. e. variation of the electric permit- 
tivity of the medium resulting from the effect of an electric polarizing field: Ae’. 

The first experimental investigation was carried out by Herweg (1920), Kautzsch 
(1928) and others. Dipolar liquids, according to their particular properties, exhibit 
either „normal“ electric saturation i. e. a decrease in the electric permittivity € under 
the effect of a strong electric field, Je < 0, or „inverse“ electric saturation, de® > 0, 
an effect first observed in nitrobenzene by one of the present authors (see A. Piekara 
and B. Piekara, 1936). This latter effect has but recently been proved to exist in some 
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other dipolar liquids, owing to the considerable experimental difficulties involved 
(see A. Piekara and A. Chetkowski, 1956; A. Piekara, A. Chelkowski and S. Kielich, 
1957; A. Chelkowski, 1958). Herweg (1920) was the first correctly to compute the 
electric saturation in dipolar gases from Debye's theory (1912). The theory was deve- 
loped to account for anisotropic polar molecules by Debye (1925) and van Vleck 
(1932, 1937), and, subsequently, for molecules exhibiting non-linear deformation. 
by A. Piekara (1935, 1937). 'The theory accounting for polar liquids is due to Debye 
(1935), A. Piekara (1939 a, b, 1947, 1950), Peterlin and Stuart (1939), Anselm (1944), 
Frenkel (1946), Booth (1951), O’Dwyer (1951) Buckingham (1956 a), Schellman (1957) 
and the present authors (see A. Piekara and S. Kielich, 1957, 1958; S. Kielich, 1958). 
However, the effect of inverse electric saturation found in certain polar liquids is 
dealt with in the theories of the present authors (1939 a, b, 1950, 1957, 1958) and 
a mention thereof is made by Schellman (1957). 

2. Electric saturation in a magnetic field, i. e. variation of the electric permittivity 
e resulting from the effect of a magnetic polarizing field: Ae”. 

This effect has not been observed experimentally in liquids; however, it has been 
shown to exist in liquid crystals by Jezewski (1924, 1926, 1929) and Kast (1924, 1927). 
Experiments aimed at detecting this effect in pure liquids by one of the present 
authors and M. Scherer (1936) only settled an upper limit for Ae”. The theory for 
a paramagnetic gas is due to van Vleck (1932), and, for a diamagnetic gas of non- 
linearly deformable molecules — to A. Piekara (1935). Recently, the theory of this 
effect accounting for diamagnetic liquids was given by Buckingham (1957) and the 
present authors (A: Piekara and S. Kielich, 1957, 1958). _ 

3. Electric saturation in an optical field, i. e. variation of the electric permittivity 
e resulting from the effect of an optical polarizing field: Ae’. 

As yet, this effect has not been experimentally detected. However, it is proved 
theoretically (see 8 6, conclusion), that using modern pulse techniques, the effect 
could be made accessible to detection. 


II. Effects of magnetic saturation 


4. Magnetic saturation in an electric field, i. e. variation of the magnetic permea- 
bility u resulting from the effect of an electric polarizing field: Ap’. 

The theory of this effect, for the case of a paramagnetic gas, was dealt with by 
van Vleck (1932), whereas that of diamagnetic liquids was the object of studies by 
Buckingham (1957). 

5. Magnetic saturation in a magnetic field, i. e. variation of the permeability A 
brought about by a polarizing magnetic field: Au”, and 

6. Magnetic saturation in an optical field, i. e. variation of the permeability u 
resulting from the effect of an optica] polarizing field: Ay’. | 

Neither of these effects has as yet been observed in dielectric and diamagnetic i 
media, because of their extreme smallness. 


a Da, Aen, en EBENEN 9 
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III. Effects of optical saturation 


7. Optical saturation in an electric field, i. e. variation of the optical refractive 
index n resulting from the effect of an electric polarizing field: An’. 

This effect is usually investigated in practice under the form of electric birefrin- 
gence or Kerr's electrooptical effect (1875). The quantity An? itself, which is known 
by the French term ,retard absolu*, is also the object of experimental investigation. 
Theories of Kerr's effect for gases are due to Voigt (1901), Langevin (1910), Born (1918, 
1933), Gans (1921), Servant (1943, 1950) and Buckingham and Pople (1955, 1956), 
whilst theories dealing with the liquid state have been given by Müller (1936), Friedrich 
(1937), Piekara (1939, 1950), Peterlin and Stuart (1939), Anselm (1947), Buckingham 
(1955), Buckingham and Raab (1957), and the present authors (1957, 1958). 

8. Optical saturation in a magnetic field, i. e. variation of the optical refractive 
index n resulting from the effect of a magnetic polarizing field: An”. 

It is the relative effect, i. e. the magnetic birefringence, as detected by Cotton 
and Mouton (1907), that is experimentally investigated. The theory of the Cotton- 
Mouton effect in liquids is due to Raman and Krishnan (1927), Piekara (1939, 1950), 
Peterlin and Stuart (1939), Snelman (1949), and Buckingham and Pople (1956). 

9. Optical saturation in an optical field, i. e. variation of the optical refractive 
index n resulting from the effect of an optical field: An°. 

Recently, Buckingham (1956 b) pointed to the possibility of optical birefringence 
being induced in an isotropic medium by a light beam of very great intensity, and 
gave a theory of the effect. 'The latter, as yet, has not been observed in practice. 

It is worth noting that, of the nine effects enumerated, only three, namely, the 
first, seventh, and eighth have been experimentally detected and investigated. The 
others still remain to be discovered. | 

The present authors have thus far published theories of five of these effects, 
namely a theory of electric saturation in an electric, magnetic and optical field and - 
of the effects of Kerr and of Cotton Mouton, respectively (A. Piekara and S. Kielich, 
1957, 1958 a, b). The present paper brings a unified, statistical-molecular theory of 
all the nine effects for dielectric and diamagnetic isotropic media. The medium is, 
quite generally, assumed to be a condensed one, composed of anisotropic polar mole- 
cules of arbitrary symmetry, undergoing non-linear deformation in an external 
polarizing field. The interaction of the molecules of the medium is determined by 
intermolecular forces of an arbitrary nature. In the first place, the so-called molar 
constants of the various effects will be derived. In turn, the constants will be related 
to the variations of the quantities e, u and n given by eqs. (1.2) — (1.4). Finally, 
the variations of e, u and n in nitrobenzene in the six effects as yet not detected are 
numerically assessed. WË en í ; 

The theory of all these effects, for anisotropic and polar molecules having the 
axial symmetry, was presented at the 7-th Colloque Ampére in Paris, 1958 
(Piekara and Kielich, 1958 c). M. 
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2. The molar constants of electric saturation 


Electric saturation in an electric field. Consider an isotropic medium of volume 


V containing N identical molecules. The spacial distribution of the molecules with ' 


respect to one another is given by the set of configurational variables t = T (r, œ), 
with r = (rj, T5, ... Py} denoting the positional variables and w = {w}, Wg, ... Oy} — 
the orientational variables of the N molecules of the medium. The latter possess 
a permanent electric dipole moment and are electrically, magnetically and optically 
anisotropic. 

Assume the medium under consideration is acted on simultaneously by an 
external electric maesuring field E" and an electric polarizing field E. In general, 
the vectors subtend the angle Q,, given experimentally. The fields applied result in 
electric polarisation of the medium as a whole, residing in electric polarisation of the 
component molecules. As the medium is assumed to be a condensed one, the distances 
separating the molecules are small, and the latter, having undergone polarisation in 
the external field, interact by their total (i. e. permanent + induced) electric moments, 
giving rise to an additional internal field in the dielectric. This additional internal 
field superposes itself upon the external field to form the so called local field existing 
within each molecule of the dielectric. 

Immersed in the local field, the molecules of the medium not only undergo linear 
polarisation (determined, in general, by a second order tensor), but are also subject 
to non-linear deformation (given by tensors of the third and fourth orders) consisting 
in deformation of the electronic shells and in some bending of the atomic bonds. 

Hence, the total potential energy possessed by the arbitrary p-th molecule of 
the dielectric in the presence of the measuring Æ’ and polarising E external electric 
fields is determined, in tensor notation, as follows: 


u(t, E", E? = u (1)® — [ur ptus Ee pope E n 


À , 
DRAG Ek «| eg 


: (a dap Eg tz 


1 [= e Ivan, , , 
Alp pang: FP PEP +.) PEP... e 


wherein u (7)? denotes the potential energy of the p-th molecule in the absence of 
the external fields (E? = E = 0), uw — the permanent electric dipole moment of 
the p-th molecule within the dielectric, 0 — the electric polarizability tensor of 
the p-th molecule, and bes), dp — tensors determining the non-linear electro-electri- 
cal deformation of the p-th molecule of the medium, to be termed the tensors of 
electro-electrical deformability of the molecule. The above tensors are symmetrical 
in all summation indices i, j, k, I, which assume the values 1, 2, 3. F’(P and F® denote 
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the local electric measuring and polarizing fields, respectively, which act upon the 
p-th molecule of the dielectric. 


The total electric dipole moment of the p-th molecule is given by the relation 


À 2 (p) 
m; (v, E^, Ey» = z GE E) (2.2) 
dF, 
which, with respect to eq. (2.1), may be rewritten as follows: 
€ H 1 ee 
m: (t, E, E? = uP + af? FP + El b FP FS) + 
e l ee 2 
x (as Uo JF FR CC A FEA oo (2.3) 


The electric local fields FX? and F@ in eqs. (2.1) and (2.3) acting on the p-th molecule 
of the dielectric consist E the | electric fields E’ and E, respectively, and of 
the internal field arising from the interaction of the molecules of the system with 
one another and with the external electric field. The problem, if restricted to dipole 
interactions, may be stated as follows 


N 
FO = E; — > D mf, FO = gt > Don mee), (2.4) 


g=1 q=1 


with mee) denoting the electric dipole moment of the q-th molecule given ‘by eq. (2.3), 
and 


1 
TH? = ET (2 dor — 9Tpai Tpai) (2.4a) 
pa 


— the dipole interaction tensor, non-vanishing for p q and vanishing for peg, 
Ty, = T, — Tp — the vector connecting the centers of the. p-th and g-th molecules 


(r, fe Yr, ve their radii-vectors) and A, — the unit tensor, x 
Iso, 
VET 5 for i Aj. « (2k) 


The next fundamental point of the theory is to calculate the component in the 
direction of the external electric measuring field, Æ’, of the vector of the electric 
polarisation P“ arising in the dielectric from the effect of fields E’ and E, which will 
be denoted by P%,. Applying statistical mechanics to calculate Pf, and accounting 
for the fact that the dielectric consists of N molecules of one kind, one may write 


D me (2.5) 


wherein 


$ mt, = mt (v, E’, E)? o? (2.5a) 
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is the component of the electric moment of the molecule in the direction of the measu- 
ring field E^, whereas m; (t, E^, E) denotes the electric moment of the p-th molecule 
given by eq. (2.3), and the a;? are the directional cosines of the angles subtended 
by the axes of the molecular reference system attached to the p-th molecule and the 
direction of the measuring field E. 

The symbol <m%> pg in eq. (2.5) serves to indicate the mean statistical value 
of the quantity m4, , computed when the dielectric immersed in the external fields E^ 
and E is in thermodynamical equilibrium. By classical statistical mechanics, for 
a condensed medium, (m) pg is defined as follows: 

N 
X we, E’, EXP) 
p=1 
TL mte E. EY? me H ` de des. din 


(mpg, = N = 26) 
2 un E, E? 


ff. fe erh GES 


wherein k denotes Boltzmann's constant, T — the Kelvin temperature of the system, 
dv, — the element of the configurational space of the p-th molecule of the system; 


N D 

u(t, E EY? with u (t, E E)? from eq. (2.1) denotes the total potential energy 
p=1 d 

of the system in the presence of the fields E’ and E. In the absence of the external 


N 
fields, >) u (v, 0) = Uy is the potential energy of interaction of the N molecules 
p=1 
of the medium. In general, Uy consists of the Lennard-Jones potential, the London 


dispersive forces, the energy of the dipolar (Keesom) and inductive (Debye-Falken- 
hagen) interaction of the molecules, and of multipole molecular interaction of various 
other kinds. 

The calculation of (mf. p g» as defined by eq. (2.6) together with eqs. (2.1) 
and (2.3), wherein the local fields are given by eq. (2.4), is generally a rather involved 
problem, and, in fact, can be solved for the linear* effects of dielectric polarisation 
and refraction (see Kirkwood (1936), Fuller-Brown (1950), De Boer, Maasen and 
Seldam (1953), Harris (1955), Buckingham and Pople (1955), Mandel and Mazur 
(1958)). In computing the non-linear effects, it is necessary to introduce simplifica- 
tions; thus, in place of eq. (2.4), the following relations will be assumed: 


H ? H 
FO Kat, FO = Fal, Re Ch) 


wherein F? and F, denote the mean local electric fields not depending directly on the 
configuration of the molecules, and o. a(? — the cosines of the angles subtended 
by the direction of these fields and the axes of the reference system attached to the 
p-th molecule of the dielectric. 

Substitution of eqs. (2.1) and (2.3) together with (2.7) in (2.6), expansion in po- 
wers of the fields F? and F, and averageing over all possible positions of the p-th 


erg wee Or ww 
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molecule with respect to the direction of the measuring field yields (see Appendix): 


e NE SON 
(MEDEE = 7 Län, I 2 os | Fit 


(pos gie o 


SON EI 
+ een. fa jj Ont + t p QE aie? on > 


90 


IR 
eg > (BoR? wf? — ôy 84) + 


JEN 7 
dee r 
ET > (Bog? of? — "c jo?) S Ji 


Tm w m A | 7 Cafe ol ence QI of OD sig )| pri. (28) 


wherein the wee are the cosines of the angles between the i-th axis of the molecular 
system of reference attached to the p-th molecule and the j-th axis of that attached 
to the q-th, H. — the angle between the electric measuring field vector E" and the 
vector E of the polarizing field. The mean values denoted by the symbol < > without : 
subscripts are those computed in the absence of the external fields (E = E = 0). 
Thus 


(2.9) 


is the mean statistical value of the function ® (t) in the absence of the external fields. 

In eq. (2.8), the directional. cosines we … of fixing the orientation of the 
axes of different molecular reference systems with respect to one another still remain 
to be statistically averaged, account being taken of the various types of molecular 
interaction implicit in the potential energy Uy. This may be effected by the definition 
of eq. (2.9) for special cases if the symmetry and structure of the molecules making 
up the dielectric are known. 

With the statistical mean <m%,>gz, it is possible to compute the quantity 
termed *molar constant of electric saturation in a polarizing electric field" and defined 


as follows: 
Pu EE: Na GES 
Su= 3 e Le Pav iz 


N, being Avogadro's number. 
The definition (2.10) of the molar constant of electro-electric saturation S% 


ell 
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together with eq. (2.8) yields 


27 Na 
135 


SM = 


e 124; 
(1 + 2 cos? Qe) [dis + ET bu / e oe? on > a 


ee du / ox (Bah off? — Ar ôn) > =F 


4aj; Me pi / Q? (f? — A (qr) px 
d "mg OS (3wik T e + 


ij d FE 4) (rs) 5 
+ Hii he b mu of? — OR Or >)! (2.11) 


The molar constant is, in general, given correctly by eq. (2.11) for condensed 
dielectrics (condensed p and liquids) composed of polar molecules of one kind 
possessing symmetry of an arbitrary type and non-linearly deformable in an external 
electric field. In the case of a gas, the potential energy of molecular interaction may 
be neglected (Uy, = 0); the statistical mean values of eq. (2.9) then reduce to isotropic 
averages, and the following molar constant „„Sy, is obtained for a gas: 


ee 2 N, ee 
gas9 M = EH (1 + 2 cos? Qee) {Sein + 77 ETT kT Hi big + 


2 O ud 4 o € e 2 3 3» T6 
+ TF (Bai; aij — ai; a) + js pa Cais Mi M; — Gi u) — js T3 Mi A | (2.12) 


The discussion of the equations derived above will begin with that of eq. (2.12) 
for gases. The first term therein, 3 dëi, is temperature-independent and accounts 
for a purely deformational effect consisting in immediate deformation of the molecules 
of the dielectric by the external polarizing field E. This effect may be derived directly 
from eqs. (2.3), (2.5a) and (2.7) by isotropic averageing only. The second term. 
12 u, bi; (TIET, which is a deformational-dipolar one, is temperature-dependent and 
accounts for the deformational-orientational effect originally computed for gases by 
one of the present authors (A. Piekara, 1935, 1937 a b). The remaining erter 
-dependent terms, namely: the purely anisotropic term 2(3 aj aj; — aj, aj) (ET) ', 
the e -dipolar term 4(3 af; u; Mj — %; HP)(ET)-? and the purely dipolar term 
2 u; u; (kT) $, determine the effect of molecular orientation originally computed for 
gases by Debye (1925). This effect had already been introduced into the theory of 
the phenomenon of Kerr and of Cotton-Mouton by Langevin (1910) for molecules 
exhibiting anisotropy only, and by Born (1918) for anisotropic and polar molecules. 
The effect of molecular orientation consists in the following: the molecules of the 
dielectric tend to an orientation in which the axis of their permanent dipole moment 
and that of their greatest polarizability are parallel to the direction of the electric 
field applied, an effect that is counteracted by the unordered, thermal motion of the 
molecules. This results in a state of thermodynamical-statistical equilibrium in which 
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the molecules attain a degree of ordering, the privileged direction being that of the 
external electric polarizing field. 

It will be seen from eq. (2.11) for dielectric liquids that the above effects, with 
the exception of the purely deformational effect, are intimately related with the mole- 
cular interaction existing throughout the dielectric. The latter, in general, is deter- 
mined by the mean statistical values of various functions of the directional cosines 
e go of the angles between the axes of the molecular systems of reference attached 
to the p-th and to the g-th molecule of the medium. The deformational effects are 
generally insignificant and may be neglected in considering polar liquids. In weakly 
polar liquids, whose polarizability is strongly anistropic, the purely anisotropic term, 
and especially the anisotropic-dipolar term, become important. On the other hand, 
in strongly polar liquids, it is the purely dipolar term and the characteristic effect 
of molecular interaction related thereto that play the decisive part; if the energy 
of molecular interaction is sufficiently great, this effect leads to the change in the 
sign of Jet, as experimentally observed. No such change in the sign of A £^ is predicted 
by eq. (2.12), which holds for gases, and wherein the molecular interaction is not 
accounted for (see A. Piekara, 1939, 1950; A. Piekara and S. Kielich, 1957, 1958). 
The discussion of eqs. (2.11) and (2.12) with respect to the symmetry of the molecule 
will be continued in § 5. 

Electric saturation in a magnetic field. Consider a system subjected simul- 
taneously to the effect of an electric measuring field E’ and to that of a magnetic 
polarizing field H. 

The potential energy of the p-th molecule of the Tem medium is now given by 


| 1 
u(t, E. Hy? = u (v)? — = ond po ro CS 


) ) (p) 
= Lë Ar ; me Fe Fink + = po 


1 9 d Ed 7 
= (as + ; dri FQ FU) + 2 PO E e (2.13) 
wherein a?" is the magnetic polarizability tensor of the p-th molecule, bem?) and comp) — 
the tensors of its magneto-electric deformability, and F5 (? _the local polarizing magne- 
tic field acting on the p-th molecule of the system. The tensor ap? j is symmetrical 
in the indices i, j, the tensor AHH! — in j, k, and the tensor dT — in the separate 
pairs of indices i, j and k, l. 

By eq. (2.13), the following expression is obtained for the electric moment of 


the p-th molecule in the presence of the fields E’ and H: 


Lu 
mj (v, E’, Ey? = pi? + KT Fm Fine + 
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As previously, the statistical average (m? pH 18 computed from eqs. (2.13) and 
(2.14), and a definition quite similar to that given by eqs. (2.10) yields the molar 
constant of electric saturation in a magnetic field (molar constant of magneto-electric 


saturation) as follows: 
em 2n N, em E > "al (pa) 
SM = i b (45 + i Lei £ On] SE: 
é 2 
4 (3.cos? Qen — 1) E des E 2 BE 2 Bog? An — of? GE e 
per ak m M (Gef? of” — 6 bu) N pr 


genna E (30; (pa) we? — ; orf”) >|) » (2.15) 


wherein Q,, denotes the angle subtended by the vector of the electric measuring 
field E’ and that of the polarizing magnetic field H. 
For a gas, this yields 


em 2x NA 
TU M = 5135 E (cii + L4 = Ki bij) Gs 


em em 2 em em 
+ (3 cos? Qem — 1) E Ci Chi t IT (3 ui bii; — ui biij) + 
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Electric saturation in an optical field. Assume the measuring field to be, as in 
the foregoing instances, a slowly varying electric field Æ’, while the polarizing field 
is a rapidly varying electric field (whose frequency exceeds that of Debye disper- 
sion), i. e. an optical field (excluding, however, ‘frequencies within the range of 
optical absorption). The latter will be assumed to be sine variable: 6 = &, cos 27 vt, 
with d. denoting the amplitude of the vector &, v — the frequency, and t — the 
time. 

As the period v^! of the a des of the rires field is so short that the 
molecules are unable to keep pace with its variations, the “effective potential energy“ 
and “effective electric moment“ of a molecule of the medium are given by the mean 
values of (2.1) and (2.3) over an entire period of the oscillations of the polarizing elec- 
tric field applied. The definition of the time mean value of a function of the time: 

1 


EIERE | (2.17) 
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yields 


A EE NT oy. (2.174) 


and eqs. (2.1) and (2.3), on averageing over the period with respect to the rapidly 
varying field, take the form: 


p SENA H Km 2 2e 4 
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and 
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wherein a? is the tensor of the optical polarizability of the p-th molecule, beep? and 
ae — the tensor of its optico-electrical deformability, and F% @) — the ON local 
polarizing field acting upon the p-th molecule of the system. For an optically inactive 
substance, the tensor Sach is symmetrical in the indices 7, j, the tensor bip — in the 
indices j, k, and ey — in the pairs of indices i, j and k, l. The SC mean value 
of the DIR of the electric moment (2.19) in the direction of the measuring 
field E’ is obtained from eqs. (2.18) and (2.19), using a method explained in the 
Appendix, in the form: 


NM il 
im) pre La = (ao j D a Hj a2 of) H+ 
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with Q, denoting the sep subtended by the vectors Æ’ and &. 
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The molar constant of electric saturation in an optical field is defined by 


Ag NA 2) PA éi acu o «c 
SM ELM <mp)p,e — |- (ME) B, : (2.21) 
“=y DE EE, Gx EE zb? 


By (2.20), this yields 
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For a gas, the foregoing equation reduces to 
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3. Molar constants of magnetic saturation 


Magnetic saturation in an electric field. When a diamagnetic medium is 
subjected to the effect of a measuring magnetic field H’ and of a strong electric pola- 
rizing field E, the potential energy assumed by the p-th molecule of the medium 
immersed in these fields is 


en PEPPER + «Jas 


l m me à x 
= iw (b) ES bee FË = E Ze FË F DL pe JE: FP Iu E (3.1) 


wherein bac?) and ony) are the tensors of the electro-magnetic deformability of the p-th 
molecule, ana I wi denotes the magnetic local measuring field acting upon the p-th 
molecule of the gege? The tensor 57? is symmetrical in the indices i, j, and cmp) — 

in the pairs of indices i, j and k, i. 
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By (3.1), the magnetic moment of the p-th molecule of the diamagnetic medium is 


m; (t, W, Ey» = (az MONTE prev) F® + 


n me | ? 
To FA Fd + … | LEE Ne (3.2) 


By eqs. (3.1) and (3.2) (see Appendix), for a diamagnetic liquid, the molar con- 
stant of magnetic saturation in an electric field is 
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with 2,,, denoting the angle subtended by the vectors H^ and E. 
For a diamagnetic gas, eq. (3.3) reduces to 


me 27 Na me Dicet 
me me 2 me me 
+ (3 cos? Li — 1) EZ — tij H yp bii Mj — bi; uj) + 


1 m m 
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Magnetic saturation in a magnetic field. By analogy, for a diamagnetic liquid, 
the molar constant of magnetic saturation in a magnetic field is obtained in the form 


‘ 


mm 2n Na ijj 
SM = 3&8 (1 + 2c0s?9,,) 13c55; + 
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wherein o is the tensor, ET ene in all indices i, J, k, of the Miche NEIN 
deformability of the molecule, and Q,, is the angle subtended by the vectors H’ 
and H. 


For a diamagnetic gas, eq. (3.5) reduces to 


27 N, ke 
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Magnetic saturation in an optical field. Similarly, the molar constant of magnetic 
saturation in an optical field, for a diamagnetic liquid, is obtained in the form 


SH 2 kb + (3 cos? Qmo — 1) E: TI ME 


m 


dj ak (3.7 
+ Hak CY Golf? off? — yd) > |) (3.7) 
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wherein c7 is the tensor symmetrical in the pairs of indices i, j and k, l, of the optico- 
-magnetic deformability of the molecule, and Q, is the angle subtended by the vectors 
H’ and c. 


For a diamagnetic gas, eq. (3.7) reduces to 


o 27 N, mo mo mo 
cM = Se z + (8 cos? Qmo — 1) Ez — dij T 
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4. Molar constants of optical saturation 


Optical saturation in an electric field. Assume the medium immersed in a cons- 
tant polarizing electric field E is illuminated by a light beam whose rapidly varying 
electric field is the measuring field ©. The vectors E and €” subtend the angle C= 

In the presence of the fields Z^ and E, the potential energy of the p-th molecule 
of the medium is 


ult, €, E? = u (cy? — (ure $2 af FP FD + J 3 


1 | [^] oe! d oe H a 
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wherein 5? and ct are the tensors of the electro-optical deformability of the p-th 
-molecule (the former being symmetrical in the indices i, j, and the latter in the pairs 
of indices i, j and k, l), and F’® is the local optical measuring field acting upon the 
p-th molecule of the medium. 

The tensor of the optical polarizability of the p-th molecule in the presence of 


the fields €’ and E is 


d? u(t, di, Ey? 
o(p) E en ee 
By (4.1), this yields 
d 
af? (c, E) = ad? + Ut FQ H ED FE EP us (4.3) 
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wherein a2 is the tensor of the optical polarizability of the p-th molecule in the 
absence of the electric polarizing field E, and the tensors beet) and ony determine the 
direct effect of the polarizing field E on the optical polarizability of the molecule. 

The mean statistical value of the optical polarizability of the molecule in the 
presence of the electric polarizing field E is defined as follows: 


o(p) (p) yp. 

d a; (t, E) yi p 

Cop = JÍ- N : >= (4.4) 
2 ues EX?) 


Tak EE, STATIN 


here, the y® are the directional cosines of the angles subtended by the axes of the 
molecular system of reference attached to the p-th molecule of the medium and the 
direction of the optical measuring field €”, and u (v, E)®- denotes the part of the 
potential energy (4.1) of the p-th molecule dependent only on the electric polarizing 
field E. 

By (4.4), and applying the method explained in the Appendix 
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The molar constant of the optical saturation in an electric field is defined as 
follows: 


Si E 


sën Ka) — (00-9. (46) 


yielding, by (4.5), for polar liquids 
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De aki TÈ (3 oof? op? m on) + 
4- Se z (Bog? fe m of) >|} (4.7) 
and, for a polar gas, 
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Optical saturation in a magnetic field. By analogy, the molar constant of the 


optical saturation in a mägnetic field (magneto-optical saturation) is obtained for 
a diamagnetic liquid: 
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Pu Pisa be a mo) E j— hg + 
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wherein c% is the tensor symmetrical in the pairs of indices i, j and k, l, of the magneto- 
-optical saturation, and Q,,, denotes the angle subtended by the vectors & and H. 
lor a diamagnetic Pu the foregoing equation reduces to 
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Optical saturation in an optical field. Finally, the case will be considered when 
not only the measuring field €’, but the polarizing field & = £, cos 2 zt v t, too, is 
an optical field. By averageing the rapidly variable polarizing field, eqs. (4.1) and 
(4.3) now yield 
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o(p) 


wherein age) and aj?" denote the tensors of the optical polarizability of the molecule 


as SL to the ten polarizing and measuring fields, respectively, and cf?) — the 
tensor of the optico-optical deformability of the p-th molecule of the medium; the 
latter is, in general, symmetrical in the first and second pair of indices separately. 

By the method discussed in the Appendix, eqs. (4.11) and (4.12) yield the mean 
statistical value of the optical polarizability of the molecule in the presence of the 
optical polarizing field: 
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wherein ©, denotes the angle subtended by the vectors &7 and €. 

By (4.13), the molar constant of the optical saturation in an optical field, defined 
in analogy to (4.6) is, for a liquid 
DE Na | 
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and, for a gas 
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If the frequency of the measuring field vector €’ and that of the I Es 
vector & are equal or differ but insignificantly, then a; = az; and chp = cs the 
tensor cu being symmetrical in all the indices, and eqs. (4. 14) may be written in 
the form 


2 00 
SM = d fs (1 -#2 cos? 2,0) Ga = 
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Optical birefringence in an electric, magnetic or optical field. It is the molar 
constants of the optical birefringence that are usually computed in optical phenomena, 
as the molar constants of the optical saturation determine the absolute value of the 
effect arising in anisotropic medium subjected to a polarizing field (“retard absolu") 
Thus, e. g., the molar constant of the optical birefringence in an electric field, or the 
so-called molar Kerr constant, is defined as 


oe Aer N, 
Bm = = Fa (25 geg Se (a?) p, ve = 90°}: (4.18) 
or, by definition (4.6) 
By = SM, 20e = 0 SM, Qoe = 90° (4. 19) 


wherein (X^ E Ga = 0° and Sao = 0° are the values of (4.5) and (4.7) for Q e = 0° 
i.e. when the vectors & and E are parallel, while <œ°> go,, = 90° and S5, ace = 90° 
represent the respective values for Q,, = 90°. i. e. when the vectors & and E are 
perpendicular. 

The remaining molar constants of the optical birefringence are defined similarly: 
B?" is the molar constant of the optical birefringence in a magnetic field, or Cotton- 
Mouton’s molar constant, and B% is that of the optical birefringence in an optical 
field, or Buckingham's molar constant (in the paper by Buckingham, the constant 
B% is defined differently and hence differs by a factor of 1/,, see Buckingham 1956 b). 

By the foregoing definitions, eqs. (4.7), (4.9) and (4.14) yield for condensed 


media 
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For a gaseous medium, the molar constants of the optical birefringence reduce to 
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om 27 Na om om l o m o m 

gas PM = 45 a — Ch T LT (3aijaij — ay 4. (4.24) 
00 27 N oo oo E? o o o o0 

SE E LL i — do + pp Dao — af à). (4.25) 


Kerr's molar constant given for a gas by eq. (4.23) consists, in general, of four 
terms. The first of these, 3055 — 55,5 


effect (1901; 1908), and is temperature — independent. As proved by Langevin 
(1910), this effect is extremely insignificant in most substances, as compared with 


accounts for Voigt's purely deformational 


the remaining terms. It is of account in monatomic gases and in substances whose 
molecules possess the spherical symmetry (see Buckingham and Pople, 1955, 1956). 
The term 2(3 bi}; u; — bj ; uj (kT)-* accounts for Born's deformational-orientational 
effect (1933). The remaining two, namely, the purely anisotropic term (3 af aj, — 
— aj, a7) (kT) * and the anisotropic-dipolar term (3 aj, Ae uj — a, WET)? ES 
for the Langevin-Born effect of molecular orientation. In liquids (eq. (4.20)) , the 
latter two, i. e. the deformational-orientational and purely orientational effects are 
influenced by molecular interaction. The term 3 €; — cÿ”,; in the Cotton Mon 
molar constant accounts for Voigt's deformational effect, whereas the term (3 a7; a7; — 
— az, aj (T): 1 determines Langevin's orientational effect (1910). It will be seen Mm. 
(4.21) that, in liquids, this latter effect is intimately connected with molecular inter- 
action. The presence in eq. (4.25) of the anisotropic term (3 ag. ai; — üt a) (kT)? 
points to the fact that the polarizing optical field also tends to order the aniso- 
tropic molecules, thus causing the isotropic medium to become optically anisotropic, 
with its optical axis parallel to the oscillating electric vector of the optical polarizing 
field. In a condensed medium, by eq. (4.22), this effect of molecular orientation is 
modified by molecular interaction. Moreover, by eqs. (4.22) and (4.25) it will be seen 


that in the present case, too, there is a purely deformational effect described by the 


term 3 ce UN eo jj an effect that is independent of the temperature and of molecular 


interaction (see Buckingham, 1956 b). 


5. The molar constant of the nine effects of saturation, computed for some special - 
cases 


In 88 2, 3 and 4, general expressions for the molar constants of the nine satu- 
ration effects originally enumerated have been derived without any. assumptions 
whatsoever as to the symmetry of the molecules and the nature of the intermolecular 
forces. Now some special cases will be considered. 

Spherical molecules. lf the molecules of an isotropic medium possess the sphe- 


-rical symmetry, the. tensors of odd order vanish: . 


4 u; = 0, OF, = b, =... = 9 - 
whereas those of even order reduce to isotropic tensors given by the following relation - 
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ships: 


1 
ijkl = -7 © (Ôi; Ont + Or Öjı + Su Gul, 


Cu = À dj du + y (Òir dr + dr Ôr;)- (5.1) 


To simplify notation, the non-vanishing components of the last tensor, which is iso- 
tropic of order 4, will be written as follows (see Buckingham and Pople, 1956): 


Cm = X + 2y = Cy» for y — 7 = & =, 


Cj ap X = CL > for. iJ) k=l 
it Y pads : c 

Se dek? (cd) fori=h, j=l. (5.2) 
All tensors, do, Coon and Gu, as well as the corresponding scalar values a, c, cj and 
c, should be provided with upper indices explaining the nature of the measuring and 
polarizing fields, respectively, e. g. aij» a^, C kI > dtr, c't, etc. The scalar coefficients o 
and c, (with suitable upper indices) denote the electric, magnetic or optical deforma- 
bility of the molecule in the direction parallel or perpendicular to that of the polarizing 

field, the nature of which is made apparent by the second upper index. 
By (5.1) and (5.2), the general expressions derived in 88 2, 3 and 4 giving the 


nine molar constants of saturation assume the following form valid both for gaseous 
and liquid media: 
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IA = 5 {ch + 2] + (3 cos? Qom — 1) (ii — (REO 
^ 27 a 
SA = (1 + 2 cos? Q,,)c%. (5.5) 


The foregoing molar constants of saturation contain only deformational terms 
independent of the temperature and of molecular interaction. This is indicative of 
the fact that in an isotropic medium composed of spherical molecules the saturation 
arises solely and exclusively from nonlinear deformation produced in the molecule 
through the direct agency of the strong external polarizing field. 

In the present case, the molar constants of-the optical birefringences (4.20) — 
— (4.22), by (5.1) and (5.2) assume the form 


oe 27 = oe 
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Thus, an isotropic medium becomes optically birefringent although its molecules 
are spherically symmetric. In this case, the medium exhibits anisotropy only because 
the strong external polarizing field induces the anisotropy directly in the molecule. 
This is precisely Voigt's deformational effect (1901, 1908), which, it will be remem- 
bered, is extremely small as compared with the orientational effect existing when 
the molecules are anisotropic or polar anisotropic. 

Polar-anisotropic molecules. For strongly polar anisotropic molecules, the small 
deformational effect may be neglected, i. e. the problem reduces to the effect of pure 
molecular orientation which now is the only one to play an important part. The nine 
molar constants derived in 88 2, 3 and 4 now reduce to the following six: 
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PAR] N | ii 3 e d 
Su SM a (3 cos? Q,, — 1) er 2 (30%? wt? — Ô; dm) 5 + 
Z 
+ Sepa K E 2 (3080 elt). dr ois (5.9) 
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The foregoing molar constants, in general, account for nine effects of molecular orienta- 


tion produced within the condensed medium by an electric, magnetic or optical 
polarizing field. 

If the polar and anisotropic molecules possess the axial symmetry, and if the 
axis of symmetry is the 3-axis, then 


u Sn ern moses pu - * d 

Du = U = 9, Ms 0, aj; = aj = aj, = 0 for ix] 

€ ee e e Oe, AE m m 0 o Z o 
di; = 45375 dan, ji = 43475 daa, 04i — 497 433- (5.13) 


Now, by the general formula of the directional cosines of the angles between 
the axes of the particular molecular systems of reference: 
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eqs. (5.7) — (5.12) yield (see Piekara and Kielich, 1958 c) 
a Aer AN ath OS AX*u 3 
Mu xs (l er 2 cos? Dee) Era Rem + WEHEN p Rg — ai Rs) D 
(5.15) 
em ` qme _ An Na ër JU" do” uà 
SM cL SM = 135 (3 cos? 0 =< 1) ra Rem + eT? Rx > (5.16) 
60! > Qe =. 4n Na Ka SEH 
Su = SM 135 (3 cos? Q,, — 1) [EE ET Rom TS GES ny). (5.17) 
mm 87 NA Jim SC 
Si" = A (1 + 2 cos? Q (5.18) 
mo om An Na 2 = 
Sm = SM = 135 (3 cos? Qmo — 1) Rem: (5.19) 
An Na a? He 


00 LE s? TR f 
SM 135 (3c Or EE kT Rem; (5.20) 


Theory of Saturation Phenomena 463 


wherein the quantities 
pe ere e gom __ „m m po o c 
4^ — ag a, A" — a$ ah, H° = ad, — Git (5.21) 


are the anisotropies of the electric, magnetic and optical polarizabilities of the medium, 
respectively. 
The quantities Roy, Rg and Rg in eqs. (5.15) — (5.20) are respectively of the form 


1 
Rem = 5 ox (3 cos? Ora) — IS? (5.22) 


> (5.23) 
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with 6%” denoting the angle subtended by the axes of symmetry of the p-th and q-th 
molecules of the medium, while the mean values symbolized by the brackets — > 
are those defined by eq. (2.9). The quantities Roy, Rx and Rg account for the effect 
of the molecular interaction of molecules possessing the axial symmetry in a medium, 
and will be termed correlation factors. 

Thus it results for the important case of a medium the molecules of which are 
axial symmetric that the effect of the molecular interaction in the nine molecular órienta- 
tional effects under consideration is accounted for completely by the three correlation 
factors Roy, Rg and Ry. This is a fact of great significance in comparing the theory 
proposed in the present paper with the experimental results (see the following § 6). 
The correlation factors Roy, Ry and Rg in the form of eqs. (5.22) — (5.24) have been 
derived and discussed for some special models of molecular interaction in earlier 
papers (see Piekara and Kielich, 1957, 1958 a, 1958 b; Kielich 1958). 

In particular, for the model of dipole pairwise interaction originally proposed 
by one of the present authors (see Piekara 1937, 1939, 1950), the correlation factors 


assume the form 
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Roy = 23 35e copo di D 
CM 7 y K y 
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wherein L = L(y) is Langevin’s function, and y = W/kT the potential energy 
of interaction of a pair of dipoles in kT units. The upper and lower signs account 
for the cases of a pair of dipoles tending to the parallel and antiparallel orientation, 


respectively. 
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It should be stresses that the foregoing model is in itself adequate for explaining 
theoretically the change in sign of the electric saturation effect produced by an external 
electric field in polar liquids, as experimentally observed (A. Piekara, B. Piekara, 
1936). Namely, in strongly polar liquids, the purely dipolar term predominates; 
then, by (5.15), 


See Ast Na 
dip Ma, 


(1 + 2 cos? Que) Rs. (5.26) 


k3 Br 
If the interacting dipoles tend to the antiparallel orientation, and for the value y > 1.33 
of the interaction energy, the correlation factor R changes its sign becoming negative 
and, hence, the molar constant of saturation (5.26) becomes positive. 

If the potential energy of the molecular interaction may be neglected, Uy = 0, 
we have 
and eqs. (5.15) — (5.20) assume a form that holds for gaseous media. 


6. Relation between the measured variations of & y and n and the molecular 
constants 


General relationships. By the first of eqs. (1.1) and by (2.5), the following general 
equation for the electric permittivity of a medium as measured by the measuring 
field E’ in the [Me of a strong electric polarizing field E is obtained: 


H 
in TE (6.1) 


wherein (mz denotes the mean statistical value of the component of the electric 
moment of a molecule of the dielectric in the direction of the field E as defined in 
general by eq. (2.6) and given, on computing, by eq. (2.8). 

Thus, for polar liquids, eqs. (6.1) and (2.8) yield 
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with a^ = + af, denoting the mean electric polarizability of a molecule in the medium 
in the absence of a polarizing field. 

In particular, the electric polarizability of a medium, measured in the absence 
of a polarizing field (E = 0), is obtained in the form 


a SEN PAC 
lie (at + 3kT S. oif s 


Eqs. (6.2) and (6.3) together yield the variation of the electric permittivity of 
the medium as arising from the effect of a strong electric polarizing field E, i. e. 
Ae’ = & — e. Restricting to small values of Ae‘ and neglecting electrostriction 
and the electro-calorie effect, by eq. (2.11) the true effect of saturation is 


eil , (6.3) 
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and Lares = denoting the molar volume (M — the molecular weight and 9 — the 


density of the dielectric), while S% is the molar constant of the electric saturation in 
an electric field, as given by the general eq. (2.11). Thus, a general expression relating 
Aet, and the molar constant S% in the case of electric saturation in an electric field 
has been derived. 

The expressions relating the variations of e, u and n given by eqs. (1.2) — (1.4) 
and the respective molar constants computed in §§ 2, 3 and 4 are derived by analogy 
for the remaining eight effects. 

In this way, neglecting the terms in the fourth and higher powers of the EE 
fields, the following set of equations accounting for all nine effects is obtained: 


A&,—S*"ENO Aen — SHY, Ae, = SC, (6.6) 
El ibt Lie e o (6.7) 
A(nft,— S% E? AnD”, = S" H. A(n), =S 22, (6.8) 


wherein specific constants of saturation S**, S@”, ... S"? have been introduced, which 
are related to the molar constants S%, ... S% computed in 882, 3 and 4 by the following 
equations: 
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here, the place of the index x should be filled in with one of the letters e, m or o. 
according to the kind of polarizing field considered. The parameters o", o"* and o°% 


—1 E r 2 
oe ( 22 di E (6.10) 


are given by 
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Nao ke AN, 

Sin AR y "3g" fo = nga 25" (6.13) 


with a” = 4 aZ and a? = } aj; denoting the mean diamagnetic and optical polariza- 
bility, respectively, of a molecule of the medium in the absence of a polarizing field; 
X is the intensity of one of the external polarizing fields, and X, — that of the local 
polarizing field arising within the dielectric through the effect of the external field X. 

Computation of theoretically predicted values of as yet undetected effects. The 
effects A(n?)°,, An), and Je, 
for some time (see Introduction). All the remaining effects given by eqs. (6.6) — 


have been detected experimentally and investigated 


(6.8) still await detection and experimental investigation. Thus, it may be useful to 
assess their respective numerical values from the present theory, using the experi- 
mental values of the constants /(n2)*,, and A(n2)”, of the known effects. The latter 
are those of Kerr and Cotton-Mouton, for which the respective constants of birefrin- 


gence K and C are defined by 


"Hin CDL np RS 
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wherein LT and n p are the values of the refractive indices for the vector of the mea- 
suring field parallel and perpendicular to that of the polarizing field. 
By eqs. (6.8), 
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alternatively, by eqs. (6.9) and (4.19), 
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with B% and B% denoting the molar constants of optical birefringence in an electric 
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or magnetic field, as given in general by eqs. (4.20) and (4.21), while o% and o” 
are the parameters determined by eq. (6.12). All parameters o°”, ... o° appearing 
in the foregoing equations may be computed theoretically from their definitions 
(6.10) - (6.12) assuming for the medium under ‘consideration the appropriate model 
of the local field existing therein. 

Assuming the model of Lorentz (1916), the local fields are of the form 


= > + 2 u +2 n? + 2 
pom E F,=——H, F,= 
à : 3 F, 3 E, (6.17) 
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If considerations be restricted to the purely orientational effect and to molecules 
possessing the axial symmetry, eqs. (6.6) — (6.8) together with (5.15) — (5.20) may 
be expressed by the constants K and C of (6.14) — (6.16). Thus, assuming the Lorentz 
model (6.17) and (6.18), the equations yielding the six effects as yet experimentally 
not detected take the form 
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where: 
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The effects defined by eqs. (6.19) — (6.24) will now be numerically evaluated 
for nitrobenzene from the experimental data at 20° C: 
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With the Lorentz model, eqs. (6.19) — (6.24) yield 
Aen. = 1.8 (3 cos? Q,,, — 1) x 102°, 


sat 


let p= 1.9 me P =) 210 


abu 
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From this set of formulas it may be concluded that, using modern equipment, 
it should be possible experimentally to detect the effects Je", and Ae?» besides the 
one, /n?,,, already predicted by Buckingham (1956). As a matter of fact, for nitro- 
benzene subjected to the effect of a magnetic field of 4 x 104 Oe parallel to the electric 
measuring field, the theory yields a value of 425, = 6 X 10-5. This represents an 
effect- accessible to measurement, though the difficulties to be overcome are very 
great considering the simultaneous by-effects (as, e. g. the rise in temperature of the 
liquid brought about by the eddy currents in the condenser plates). Such experiments 
are being performed in this Laboratory, using modern techniques. 

The second of the effects accessible to detection, As‘, amounts to as much as 
1x 10-5if £2 = 4.2x 108 erg/em?, i.e. if 2, = 19.5 kV/cm. An electromagnetic wave beam: 
whose electric vector is of such great intensity should have energy density amounting 
to e 2/8 or a luminous flux of (cn/8x) & = 109 W/em?. The latter could be obtained 
only as a short flash, by discharging a high voltage condenser battery through a gas. 


Appendix ` 


The rather involved calculations in the method proposed by the present authors 
will be exemplified by one of the consecutive steps, namely, the one leading to the 
computation of the statistical mean value of the component of the magnetic moment of 


iu 


rain 
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a molecule of the diamagnetic medium in the direction of the measuring field H^, 
in the presence of an electric polarising field E: 
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here, m7 (t, H’, EY? and u(t, H^ E)? denote the magnetic moment and the potential 
energy of the p-th molecule of the diamagnetic medium as given by ‘eqs. (3.1) and 
(3.2), respectively. The directional cosines of the angles subtended by the axes of the 
molecular system of reference attached to the p-th molecule of the medium and the 
direction of the magnetic measuring field H’ or of the electric polarizing field E are 
denoted by B® and a(?, respectively. 

Assuming, by analogy with eq. (2.7), that 


i Po p pe, p T aP, 


and on expanding in powers of the magnetic and electric fields, eq. (I) yields 
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wherein the statistical mean values symbolized e the brackets < > are those com- 
puted in the absence of external fields (H^ = E = 0), and are defined by eq. (2.9). 

In the absence of external fields, the various orientations of the p-th molecule 
with respect to the laboratory referential of the medium are all equally probable. 
Hence, in the expansión (II), it is permissible to average the directional cosines BP, 
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… af” separately and isotropically over all the possible orientations of the p- -th mole- 
cule of the medium with respect to that system of reference. To do this, the directional ` 
cosines of the q-th and 7-th molecules should first be expressed as functions of those 
of the p-th molecule by the transformation formulas 


a) = oD aD, aP = eX? a? 


in the case under consideration, when the molecular systems of reference attached 
to the various molecules are rectangular, the coefficients of the transformation: 0%”, 
wH”, have the meaning of the cosines of the angles subtended by the axes of those 
systems. Thus, e. g., w29 is the cosine of the angle subtended by the o-axis of the 
molecular systems of reference attached to the p-th molecule and the k-axis of that 

attached to the q-th. 
By the foregoing transformation formulas and by eq. (5.14). the isotropic mean 


values of the functions of the directional cosines are 


aP = pO = = pe p a = =... =Q, 
(ei porc» jl OBI 
f CONS 


r) Ce l rs 
SE BP gi af (r) (Da? = =z {2 (2 eos? Que) wif? wy E 


+ (3 cos? Qne — 1) (08 c 49 -L ef? of}, (I) 


wherein A, is the unit tensor of (2.4b), while H. is the angle subtended by the vectors 
H and E. 
With (III) and eq. (5.14), the expansion UD takes the form 
1 
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With the molar constant of magnetic saturation in an electric field defined as 
follows: 


me _ © 
td SM = -7 p s (mmp?m,E— = (mmp»m, 2) N (V) 


eq. (IV) yields eqs. (3.3) and (3.4). 
The molar constants of the remaining eight effects are computed similarly. 
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Three methods of obtaining ferroelectric systems with a small temperature coefficient of 
the permittivity TC e are presented. The first consists in connecting in parallel two fer- 
roelectric condensers whose TC e differ in sign, the second in composing heterogeneous 
systems from two or three different BaTiO, — SrTiO, solid solutions, and the third in 
composing BaTiO, — MgSnO, solid solutions. The last method presents the best advan- 
tage yielding | TC e| < 200 x 1074 deg”! and a rather high value of the permittivity. The 
TC &( T) dependence near the ferroelectric transition temperature is discussed. From the 
dielectric properties investigated, the BaTiO, — MgSnO, solid solutions reveal a beha- 
viour typical of other perovskite type ferroelectrics. 


Introduction 


The discovery of the anomalous dielectric properties of Rochelle salt by Valasek 
(1922) and the establishment of the analogy between these properties and those of 


ferromagnetic materials gave rise to a new class of dielectrics known at present as 


the ferroelectric materials. The discovery of the ferroelectric properties in potassium 
dihydrogen phosphate by Busch and Scherrer (1935) and in barium metatitanate 
by Vul and Goldman (1945) and, independently, by von Hippel and his coworkers 
(1946), were milestones in this field of investigation. Recently, Matthias and Remeika 
reported on dicalcium-strontium propionate (1957) which seems to be a representant 
of a new class of ferroelectrics. Except for the latter, ferroelectric materials can be 
classified in two groups. The first contains crystals with so-called hydrogen bonds 
of various kinds, whilst the second contains crystals with oxygen octahedra. 
Ferroelectricity results when spontaneous polarization arises within a crystal; 
this may occur as a result of dipole interaction. The acting dipoles may arise from 
off-centered protons in the H-bonds of the first group of ferroelectrics, or from off- 
-centered metal ions of the oxygen octahedra of the second group. The existence of 
spontaneous polarization in the crystal involves unusual properties which may occur 
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only in a number of crystallographic structures, as they are closely related to the 
properties of the crystal lattice. Thus, spontaneous polarization may exist only in the 
10 pyroelectric crystallographic classes possessing no symmetry centrum. The spon- 
taneous stresses involved by the spontaneous polarization cause a lowering of the 
symmetry in passing through the transition temperature from the paraelectric to the 
ferroelectric phase, determining the crystallographic class of the crystal in the Curie 
region. The theory of ferroelectric transitions proposed by Zheludyev and Shuvalov 
(1957) makes it possible to predict within which crystallographic classes this pheno- 
menom can take place. Sufficiently high temperatures counteract the spontaneous 
polarization, and cancel the dipole interaction above the Curie point. The ferroelectric 
becomes a paraelectric and the permittivity decreases, following the well-known 
Curie-Weiss law as expressed by the formula 


G 
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where £ denotes the „background“ permittivity, c, and 0, being the Curie constant 
and Curie-Weiss temperature, respectively. At the same time, 0, denotes the tempe- 
rature of the Mosotti catastrophe for the paraelectric phase. Similarly, the permittivity 
of the ferroelectric phase for the temperature range near the transition point follows 
an analogous Curie-Weiss law 
Cf 
E= D += 
Ra GC 

with a different Curie constant and Curie-Weiss temperature. 

A satisfactory theory of ferroelectricity does not exist as yet. For the most impor- 
tant representant of the oxygen-octahedra group of ferroelectrics i. e. barium metati- 
tanate having perovskite type structure, a number of theoretical treatments endeavour 
to describe the main properties. Among these, Devonshire’s thermodynamical theory 
(1949, 1951) and Mason-Matthias’s statistical theory (1948), are most useful. The 
latter theory was recently generalized by Piekara (1958). 

Devonshire’s free energy function 

iie Pet uu EE 
2 4 (OS 
wherein & and ¢ are temperature-independent coefficients, implies knowledge of the 
e (T) function as following the Curie-Weiss law for the paraelectric phase. When 
the coefficients y, é, ¢ are known, the quantities 
9A An E 
E =. and — = —— 
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may be computed tor the ferroelectric phase assuming continuity of the constants 
throughout the paraelectric and ferroelectric ranges of temperature. The discontinuity 
of the first order transition, however, does not fully justify this assumption. 
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The theory of Mason and Matthias considers the permanent dipoles y arising 
by Ti ions displacement from the center of the unit cell. The internal field is 


Dd duc Ee 


The dipole polarization in a form taking into account the difference in well energy 


(W/kTA 0) and the anisotropy of the polarizability (x, = &,) is given by Piekara 
as follows: 


noc sinh: Z 
> .coshZ a 


where 


Z = uFfkT and œ = 2 exp | a : («4 — a.) cl 


W denotes the difference between the potential energy of the wells, perpendicular 
and parallel to the c-axis, in which the Ti ions possess minimum energy. Accounting 
for unit cell deformation by spontaneous polarization P, for the free crystal and using 
the relations: 


WIET = cP2. and’ p= el + AP), 


Piekara derived a new function P,(T) revealing a first order transition, with two 
critical temperatures 0, and 0, for which Le = 0 i.e. the Mosotti catastrophe occurs. 
Thus, the latter theory presents an advantage particularly in the description of the 
ferroelectric behaviour at temperatures near the transition point, which is of great 
theoretical and practical significance. Since the highest permittivity values observed 
lie between 0, and 0, (which may be termed the Curie-Weiss region), knowledge 
of these limits is highly important. 
The temperature coefficient of the permittivity 


for the paraelectric phase is determined principally by the Curie-Weiss law, and that 
for the ferroelectric phase — by the ferroelectric Curie-Weiss law already mentioned. 
In investigating various ferroelectric materials with the intention of choosing those 
presenting high values of e with low TCe values, the general rule determining the 
Curie-Weiss law for different compositions and preparations of-the material to be 
tested should be known; however, no such rule is known as yet. Thus, the results in 
the present paper have been obtained without the indications which would have been 
given by such a theory. The present investigation was intended to answer the question 


. whether it is possible to produce ferroelectric systems of the perovskite type with 


small TC e, and to point to the principle upon which the solution of this problem 
should rest. ds m i 
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Samples 


Samples of BaTiO, and of BaTiO, — SrTiO, solid solutions with different compo- 
sitions were obtained following up numerous data on this problem (Bunting et al., 
1946, Smolenski, 1950). In turn, a technology of production of heterogeneous systems 
of BaTiO, — SrTiO, solid solutions and homogeneous systems of BaTiO, —MgSnO, 
solid solutions was elaborated. 

Attention was directed to such problems as the presence of impurities, the in- 
fluence of sample support on the process of preparation, and the grain size of the 
material. Obviously, other decisive factors such as pressure as well as the temperature 
and time of sintering should be taken into account in the synthesis. These factors 
jointly determine the behaviour of the specimens obtained. 

The importance of impurities in producing the samples from chemically pure 
reagents should not be underestimated. Although opinion is contradictory (Vul and 
Goldman, 1948, Eubank, Rogers, Schilberg and Skolnik, 1952) as to the róle played 
by impurities in ferroelectrics of the type in question, it seemed desirable to know 
the contents of the impurities in the material used. The spectrographic method was 
employed to determine the principal impurities in the reagents and in the intermediate 
and final products of reactions. The impurities found were: Si, Mg, Ca, Fe, Sn in titania 
(TiO,), Sr, Ca, Mg in barium carbonate (BaCO,) and Ba, Ca, Mg, Si, Sn in strontium 
carbonate (SrCO,). Comparison of the dielectric behaviour of samples having somewhat 
different contents of impurities did not reveal any essential differences between them. 

The influence of sample support on the properties of the specimens obtained 
is very considerable. High sintering temperature affected the reaction rate between 
the sample and the material of the plate on which it rested during the production 
process. Any reaction occuring leads to additional impurities. Platinum, alumina 
(A1,0,) tablets or plates made from the same material as that obtained in the process 
of synthesis yielded the best supports. The latter rendered the highest degree of 
purity, completely excluding the introduction of impurities during the sintering 
process. Moreover, an attempt was made to use other materials as supports; however, 
the reaction of the sample during sintering with such materials was much too rapid 
at the temperatures employed. This occurred, for instance, when using an oxidated 
nickel plate, which caused remarkable changes in the properties of the samples owing 
to diffusion of foreign atoms. Thus, the permittivity of such contaminated ferroelectric 
material decreases in certain conditions to about fifty percent of the value measured 
in the pure material. 

The grain size of the reagents is of importance for the reaction rate in solid state 
substrates. Thus, the reagents were milled to the smallest grain size in order to develop 
the largest active surface. This was performed in a porcelain mortar to a grain size 
of less then 50 microns in diameter. The grain size of the intermediate and final pro- 


ducts depends on the milling time and on the pressure and sintering conditions of - 
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the process. To avoid high porosity of the material and introducing inhomogeneity 
from large microcrystals in the ceramics, the grains should be small. On the other 
hand, it should be expected (cf. Jaccard, Kanzing, Peter, 1953), that lowering the 
grain size below a limiting value may involve the vanishing of ferroelectric proper- 
ties of the material owing to the energy requirements of the single domains of a ferro- 
electric. Little can be seid about this question, as after milling two fractions were - 
separated and the properties of only two kinds of samples, namely, one made from 
a 50—150 microns fraction and one of less than 50 microns were compared. The main 
difference observed between them consisted in a somewhat grater value of the permit- 
tivity (about 10 per cent) in the sample obtained from the material with the smaller 
grains. Such a sample also appeared to have a lower porosity, and observations made 
with a metallographic microscope showed that the grain size of the finished sample 
was then smaller as compared to those made from the 50—150 microns fraction. 
This last fact was obviously to be expected if the other parameters of the synthesis 
were the same for both fractions. 


Obtaining small porosity of the samples required large pressures. Thus, attempts 
were made to compress the material up to the very high pressure of 35000 kG per cm?. 
It was observed that by raising the pressure from 1000 to 10000 kG per cm? one 
can obtain a lowering of the porosity which leads naturally to a rise in the effective 
value of the permittivity by about 80 per cent. The pressure was exerted with the 
aid of a large oil hydrostatic press of up to 100 tons!. The milled material was compres- 
sed in a special matrix made from improved high quality tool steel. Neither slide 
materials nor lubricants were added, to avoid introducing new impurities in the synthe- 
sizing samples. There were even difficulties which required the elimination of the 
part of the sample material contacting the side of the matrix, since it became polluted 
with iron therefrom. Such and other difficulties involved caused limiting of the pres- 
sure to about 1000 kg per cm? in most cases. 


The important factors determining the behaviour of the samples produced are 
the temperature and time of the sintering process. As a general rule it may be assumed 
that, up to about 1350°C, the higher the temperature and the longer the time of sin- 
tering the higher is the permittivity obtained, this being due chiefly to the lowering 
of porosity of the samples. As shown by Jezewski and Piech (1955), repetition of the 
sintering process in the same material also leads to an increase of the permittivity. 
Nevertheless, precaution should be taken if the sintering temperature is to be raised 
above 1350°C. The danger consists in the possibility of the transition of barium 
metatitanate from the ferroelectric perovskite into the non-ferroelectric hexagonal 
structure, which takes place at about 1450°C (Trzebiatowski et al., 1954). Longer 
sintering times do not. influence the direction of the reaction; however, a limit 
was set by the thermal ‘endurance of the laboratory furnace employed. In appro- 
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priate conditions the synthesis of metatitanates generally follows the equation: 
MeCO, + TiO, > MeTiO; + CO; 


where Me stands for the bivalent ion e. g. Ba** or Sr++. This indicates that the material 
obtained must be porous owing to carbon dioxide escaping during the reaction in 
the solid state. Hence, the synthesis was carried out in two stages in order to avoid 
this disadvantage. This method allowed of a shortening of the sintering time, which 


could be limited to several hours. 


a) Stage I 

After 3 hours drying at 150?C, the constituents of the batch weighed in the 
molar ratio were thoroughly milled together for 1 hour in a porcelain mortar to a grain 
size of less than 50 microns. The mixture was then compressed without any addition 
in a steel die at 1000 kG per cm? into tablets of a diameter of 13.6 mm, the weight 
of a tablet being 2 grams. The tablets were then heated at 1300°C for 3 hours in an 
electric muffle furnace with sylite heating resistors. The desired temperature was 
attained in 1 hour and measured with a Pt-Pt/Rh thermocouple with an accuracy 
of + 20°C. The tablets rested on a barium-strontium plate supported by alumina 
tablets or platinum wires. On cooling overnight, the first stage of the preparation 
was completed. According to investigations (Trzebiatowski and coworkers 1950, 
1952, Wojciechowska, 1953) on the mechanism of the synthesis of barium and stron- 
tium metatitanate, it may be concluded that the porous product obtained in this 
way consists of BaTiO, or of a solid solution without any instable transition compounds. 

b) Stage II 

Next, the porous tablets were milled and passed through a sieve having 10000 
meshes per cm?. Samples were then made by compressing the powder once more 
at 1000 kG per cm? into tablets of a diameter of 13.6 mm and weighing 1.6 grams, 
and refiring them at 1200°C for 2 hours and at 1300°C for 3 hours. In some cases, 
moreover, higher pressures amounting to 20000 kG per cm? produced in a special 
steel die were employed. The hard ceramic samples thus obtained were suitable 
for investigations, as well as for the further preparation of heterogeneous mixtures. 
The good samples were about 7 grades hard by the Moss hardness scale and frequently 
possessed a glassy and shining breaking. 


ec) Stage III 


To prepare the heterogeneous mixtures described below, solid solutions recrystal- 


lized during the second stage of preparation were crushed and passed through sui- _ 


table sieves. Two fractions of the material powdered were u&ed. The grain size of 
the first, measured with a. microscope micrometer eyepiece, ranged from about 50 
to 150 microns, being less than 50 microns for the second. Two, and in some cases 
three different solid solutions of the same grain size were used to prepare a hetero- 
geneous system. After drying, weighing and mixing thoroughly for 1 hour, new 


tablets of the above-mentioned dimensions were obtained by compression at 1000 kG- 
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per cm? and sintering in an electric furnace during 4 hours at 1000°C and 2 hours 
at 1200°C, so that the reaction between both constituents could proceed only within 
the boundary layer of the grains. 

Homogeneous MgSnO, — BaTiO, systems were also prepared in three stages 
in order to obtain low porosity of the material in question. In this case, the third 
stage was a repetition of the second one; the sintering temperature during the last 
half hour was, however, maintained at 1380°C. 

To obtain a regular shape of the samples, both faces of the tablets were thoroughly 
ground and polished until exactly parallel. The samples were weighed and measured 
with a micrometer and subsequently coated with silver paste and heated to 700°C 
for 15 minutes to form silver electrodes on the faces. The silver paste consisted chiefly 
of silver oxide (Ag,O), thoroughly milled together with glycerol in a porcelain mortar. 
As terminals of the ferroelectric condenser obtained, thin copper wires were then 
soldered to the silver electrodes. Moreover, some of the samples were coated with 
silver electrodes by way of evaporation in vacuo. For this purpose, a special mask 
was made, enabling a layer of silver to be deposited simultaneously on 25 tablets. 
After evaporation, the tablets were investigated in a suitable holder constructed with 
a good invariable insulator. 

For preparation, reagents of the c. p. quality produced by The British d Houses, 
Ltd., Poole, England (titanium dioxide) and by The Laboratory Chemicals Factory, 
Gliwice, Poland (the other compounds being barium, strontium and magnesium 
carbonates and stannic dioxide) were used. X-ray. verification was carried out to 
eliminate TiO, reagents exhibiting the brookite or anathase modification, as their 
presence would complicate the process of synthesis by necessitating a trarisition into 
the stable rutile form. In this way it was found that the reagents used consisted enti- 
rely of the rutile modification. In the ready products of the reactions, only a perov- 
skite structure was roentgenographically observed?. This meant that the process of 
the synthesis had proceeded completely in the desired direction. 

A number of samples were produced from additionally purified reagents of the 
above-mentioned quality. The aim of this operation was to eliminate compounds 
occurring in the substrates of the synthesis, such as certain oxides, chlorides and 
sulphates, since the presence of these impurities might influence ‘the resistance and 
electric breakdown of the samples produced. Thus, it seemed worth while to purify 
the reagents by way of multiple washing with distilled and, subsequently, with redis- 
tilled water. The titanium dioxide was previously treated with hydrochlorid acid in 
order to eliminate impurities soluble therein. The water used in the elutions was 
analysed for traces of chlorides and sulphates and its resistance was measured. When 
the salt traces had ceased to be detectable and resistance was established at a constant 
level, the process of purification was considered to be completed. Measurements on 
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samples produced from such purified materials showed conclusively that their resi- 
stance rose by about one order of magnitude as compared to the value obtained with 
samples from non-purified substrates. Their electric breakdown rose by about 50 per 
cent. Thus, such an additional process of purification appears to be indicated when 
working with reagents of c. p. quality. 

Further attention was directed to the iron impurity introduced into the sample 
when compressing the material in the iron matrices, particularly on applying higher 
pressures. To avoid the polluting of samples as occurring in the part of the material 
contacting the inner walls of the matrix, a special matrix model was made layed out 
with zinc. Thus, the high friction produced by pressure on the walls of the matrix 
now involved polluting with zinc instead of iron. Zinc was chosen because of its low 
boiling temperature (near 906°C), so that the impurity thus introduced should eva- 
porate when sintering the material at high temperatures. The samples compressed 
in such a zinc matrix did not show any zinc impurity. Samples compressed in the 
zinc-coated matrix showed a somewhat higher resistance that those compressed in 
an iron matrix. For this reason, the latter samples were slightly ground preceding 
the measurements, to remove the outer layer polluted with iron. 

A comparison of the dielectric properties of samples made from materials milled 
in different mortars showed that porcelain does not introduce any noticeable impuri- 
-ties as compared to agate. Therefore a porcelain mortar was used to produce most 
of the samples examined in the present investigation. 

Lastly, the problem of the possible pollution of the samples in the sintering process 
was considered. Spectrographic checking, however, showed that no new impurities 
were introduced when sintering the samples in the manner described above. 


Measuring devices 


The capacity of the samples investigated was measured by the radiofrequency 
resonance method. The circuit consisting of a HF oscillator and a bridge serving 
as receiver is shown in Fig. 1. The tuned resonant circuit of the bridge contained 
the sample measured and the variable calibrated precision condenser with direct 
sapacity reading. This circuit was loosely coupled to a Hartley oscillator by means 
of two coils about 30 cm distant from each other, the frequency of the oscillator 
being 1 MHz. The resonance curve of the tuned circuit with and without the sample 
was determined on the bridge with a galvanometer of adjustable sensitivity. The 
first valve of the bridge acted as grid detector and amplifier (audion), while the second 
one compensated the galvanometer current to zero, using an adjustable cathode 
resistor. The measurements were carried out near the maximum on both slopes of 
the resonance curve, any variations of the capacity measured being compensated 
with adjustable calibrated precision condensers by Jezewski's method (1927, 1933), 
which makes it possible to measure the permittivity of the material in question in- 
dependently of the dielectric power factor and conductivity of the sample. The mea- 
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suring HF-tension on the ferroelectric condensers could be determined with the aid 
of the HF valve voltmeter in the receiver resonant circuit. The power factor of the 
samples could be estimated by measuring the lowering of the resonance curve maxi- 
mum. To eliminate voltage fluctuations, a valve stabilizer was used as power supply. 
Carefuly operated storage batteries supplied the filament voltage to the valves The 
apparatus was assembled compactly within a grounded shield, observing the usual 
precautions as to screening and location of the elements. The maximum sensitivity 
as well as the accuracy of the method in measuring € were of the order of 1075, the 
measuring HF-tension on the sample being 0.5 V, whereas the capacity measured 
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Fig. 1. Radiofrequency resonance method 


ranged up to 2000 pF. As only the relative changes in permittivity of the sample and 
not the absolute values were the principal object of the present investigation, the 
precautions regarding the shape of the measuring condenser as well as the inductances 
on the connections could be neglected if the measurements were carried out using 
the same connecting cables and samples of the same shape. =: 
The samples were heated in a Hoeppler ultrathermostat, or in a glycerol-thermo- 
stat especially built for the purpose with a temperature range of up to 200°C. The 
samples having no copper connections soldered on were maintained in a specially 
constructed holder which made it possible to heat and to measure up to six ferroelectric 
condensers simultaneously. Several hundred samples were prepared and investi- 
gated, and so this improvement reduced the time required for the measurements. The 
sample holder was provided with a system of thermocouples for checking the small 
differences in temperature between the precision thermometer and each of the mea- 
sured ferroelectric condensers. The thermocouples were made of thin copper and 
Eureca resistance wires (cf. Piekara and Pajak, 1952). By means of a seven-point 
rotary switch, six different samples and a checking mica condenser could be succesi- 
vely connected to the apparatus with a coaxial cable at the same measuring temperature. 
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Some of the measurements were carried out in a special vacuum thermostat, 
making it possible to cover the temperature range from — 180°C to +180°C. 

By the mixture law, the true permittivity of the samples investigated could be 
computed from the porosity and effective permittivity of the samples. Thus, it was 
necessary to determine the porosity and density of the material. The porosity was 
calculated by the following formula: 
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where d, is the pycnometric density of the powdered material, and d, — the effective 
density of the sample. The last was measured by the hydrostatic balance method 
or by measuring the dimensions of the sample with a micrometer, using samples of 
regular shapes. 

When determining the pycnometric density of the finely pulverized material, 
one should take into account the moistening of the powder with the pycnometric 
liquid. Using water as a pycnometric liquid, a considerable lowering of the measured 
values of the density was observed as compared to the known X-ray density of the 
material. Boiling the powder for various periods of time up to 7 hours in redistilled 
water wherein it was subsequently weighed in the pycnometer yielded an increase 
in the density. This probably means that air was occluded on the powder, lowering 
the measured value of the density. Finally, it was decided to carry out the measure- 
ments on pulverized materials boiled during 7 hours, since such measurements were 
found to be in good agreement with those carried out for the purpose of checking, 
using petrol or carbon tetrachlorid as pycnometric liquid. 

The dielectric power factor was estimated by measuring on the HT Schering 
bridge shown in Fig. 2. The electric breakdown was estimated for some ceramic sam- 
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Fig. 2. Schering bridge 


ples immersed in a paraffin oil bath using a direct current (up to 5 kV) and an alterna- 
ting current (up to 13 kV) power supply. This power supply was used when a biasing 


field was applied on the samples, in the way already described (Piekara and Pajak. 
1952). | | 
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Dielectric hysteresis loops were observed using the Sawyer-Tower (1935) method 
(Fig. 3). The aim of the method used is to check the existence of a ferroelectric region 
in the material investigated. This problem is of interest for heterogeneous systems 
having no distinct permittivity maximum, because the investigation of hysteresis 
loops reveals a Curie temperature, if present. 


mtm 


Fig. 3. Sawyer-Tower method 


Microscopic observations made with the aid of a metallographic microscope, 
using the micrometer eyepiece, made it possible to estimate the grain size in the ready 
sintered materials after their surface had been suitably ground and polished. As 
etching liquid, hydrochloric acid was used and etching was performed at room tem- 


perature. 


Compensation systems 


The fact that well known ferroelectric systems exhibit a strong temperature 
dependence of their permittivity does not cancel the possibility of using them as 
condenser materials. Three different methods will be presented for obtaining systems 
with a small temperature coefficient of the permittivity as compared to the values 
usually observed. The first, simple method consists in connecting in parallel two 
ferroelectric condensers whose temperature coefficients of the permittivity differ in 
sign. Considering the temperature dependence of the permittivity for the barium 
metatitanate class of ferroelectrics, one can see that they present the possibility of 
using the two parts of the e (T) curve, separated by the Curie point, for compensation 
systems. Connecting two ferroelectric condensers in parallel, it is here necessary to 
use materials possessing Curie points situated below and above the temperature 
range within which the system is expected to work. Thus, detailed investigations were 
undertaken on e (T) curves of barium-strontium metatitanate solid solutions to get 
information as to the choice of the materials needed. "d 

By the first two stages of the process of preparation already described, samples 
consisting of simple solid solutions of various concentration of SrTiO, in BaTiO, 
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Fig. 4. e and TC e (dashed line) versus temperature for BaTiO; — SrTiO, solid solutions (for concen- 
tration, see Table I) 


were obtained, exhibiting a Curie point in the interval of 0°C to 120°C. The grain 
size of the samples was estimated at about 1—5 microns on the average. 

Fig. 4 shows the temperature dependence of the permittivity and the TC e (T) 
curve derived therefrom, for one composition. Fig. 5 shows the Curie point T, as 
function of the concentration in per cent by weight. The dashed line represents the 
results of Grünicher and Jakits (1954), calculated by the present author in per cent 
by weight, and is in good agreement with the foregoing curve. It seems that the insigni- 
ficant deviations in the results obtained, consisting in a somewhat different Curie point 
for the same composition, can be explained by the difference in impurities and in 
the régimes of preparation employed. 

Fig. 6 shows the density-concentration dependence obtained from the density 
estimation made as described above, for grain sizes of the solid solutions not exceeding 
50 microns. The dashed line corresponds to the X-ray density of the solid solutions 

in question, as calculated by the present author from the lattice constants given by 
 Grünicher and Jakits (1954). The density is found to be a linear function of the concen- 
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tration (by weight). The pycnometric data are in good agreement with those calculated; 
however, the former yield values about 3.5 per cent lower. This difference is probably 
due to the absorption of air on the large surface of the powders. 
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In Table I are compared the average values of such properties as the Curie points, 
permittivity (effective and true), bulk density, porosity, electric breakdown and 
power factor (at 50 Hz, room temperature) of solid solutions tested, some of which 
have been used for compensation ferroelectic systems. The concentrations are given 
in per cent by weight; the true permittivity (i. e. on eliminating porosity) is calculated 
using Wiener's mixture law (see below), assuming the form factor u to equal 225. 

All solid solutions tested exhibited regular hysteresis loops, which, however, 
did not vanish at the Curie points measured, i. e. at the temperature of the maximum 
value of e, but persisted over a range of some centigrades more. The values of the 
spontaneous polarization, total polarization and coercive force at 50 Hz, were of the 
order of those quoted in the literature. 
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Fig. 6. Density versus SrTiO, concentration in BaTiO: — measured, - - ——- - calculated. 
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The material tested reveals a distinct dependence of the dielectric permittivity 
on the measuring field strength. Fig. 7 shows one of the e (E) responses measured on 
a 50 Hz fed HT Schering bridge at room temperature. As the Curie point is appro- 
ached, the saturation effect of the permittivity (i. e. a decrease of the permittivity 
following its maximum value) is observed at lower field strengths. Above the Curie 
point only a lowering of the permittivity values due to nonlinearity of the dipole 
polarisation occurs. 

The typical biasing field effect, consisting in an abrupt increase followed by a 
slow decrease of the dielectric permittivity, was observed in all samples. These so called 
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Fig. 7. e versus 50 Hz measuring field strength for solid solution “G“ 


“delay effects” have already been reported on earlier (Piekara and Pajak, 1952, Pajak 
et al., 1956) and are at present being investigated in detail for both short and long ob- 
servation times by Kaczmarek and Pietrzak. 

Compensation systems composed of solid solutions "B" and “M? with Curie 
points at 70°C and 0°C, respectively, have been investigated. Fig. 8 shows the capa- 
city-temperature dependence for such a “BM” condenser compensation system within 
the temperature range of 10°C to 100°C, and also the analogous responses obtained 
for the following compensation systems: “BG”, “GM” and “BGM”. The component 
condensers are prepared so as to have almost the same peak capacity of about 1000 pF. 
The same diagrams contain the TC for these systems, calculated as 
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One can see that the TC values of such compensation systems are smaller that those of 
the common simple solid solutions and yield practically useful temperature ranges of 
about 50 degrees centigrade within which ITCel < 200 x 1074 deg! (“BM” — sy- 
stem), whereas in simple solid solutions this range covers only several °C. 

It follows that the method presented makes it possible to assemble various systems 
with reduced TC values by using two or, preferably, a greater number of ferroelectric 
condensers with different Curie points. The working temperature range is easily 
modified and the TC value corrected as desired. 


Heterogeneous Systems 


The second method for obtaining a low temperature coefficient of the permitti- 
vity consisted in composing a heterogeneous system of two or three simple BaTiO, — 
SrTiO, solid solutions in the manner described above. Two solid solutions with different 
Curie temperatures but of the same size of grain are chosen and react on the grain 
surfaces in the process of sintering leading to the formation of an intermediate bound- 
ary layer. This layer consists of a solid solution having a Curie point lying between 
those of the initial solid solutions. Thus we have to deal with a material consisting of 
three or more solid phases and, moreover, a gaseous phase, owing to the porosity of 
the ceramic sample (Pajak and Piekara, 1955). The composition of the heterogeneous 
system is determined by the sintering régime and lies between the two extremal cases 
represented by a mixture of two solid phases (neglecting the gaseous phase) and 
a system of one solid phase (if the reaction has proceeded as far as to form a homo- 
geneous phase with a clearly determined Curie point). 

The results of this method are shown in Fig. 9, which brings the permittivity 
measured and the temperature coefficients of the permittivity (estimated from these 
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measurements) of heterogeneous systems composed (for concentrations, see Table I) 
as follows: 

“BG” D ep and “G” solid solutions 

“BM?” ne EE and “M” solid solutions 


“GM” 1:1 Sl and “M” solid solutions 
“BGM? — 1:31: 4B", — *G" and “M” solid solutions 


One can see that such materials reveal interesting temperature responses with no 
clear Curie points, but with high values of the dielectric pefmittivity. The optimal 
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Fig. 9. e and TC e (dashed line) versus temperature for heterogeneous systems 


temperatures for such systems range between the Curie points of the initial simple 
solid solutions. A comparison of these results with those mentioned in the previous 
section for compensating systems shows that in this case the systems obtained present 
greater advantages. The ranges within which |TCe| < 200 x 1074 deg! cover about 
100 degrees centigrade or more, and for the *GM" systems the permittivity is practi- 
cally temperature independent from 10°C to 35°C, with e = 1800. 

The results obtained may be explained by comparing the experimental responses 
with those calculated from Wiener’s (1912) or even Lichtenecker’s (1924, 1926) for- 


mulas for mixtures, namely: 


Em— 1 & —1 
———— = d vi——, or loge, = v; log &;. 
Em ar u i Ei SF H F 


which latter is applicable if the permittivities of the components do not differ consider - 
ably (Rushman and Strivens, 1947). Wiener’s formulae for a two component mixture 
with permittivities ¢, and €, reduces to : 


Em + Zës er & + Ze 
assuming the form factor u = 263. which is applicable if the second component is 
dispersed throughout the first in the form of spheres. This assumption can be made 
as a first approximation for the discussed heterogeneous mixture. Exact calculations 
cannot be made, as the form factor u for mixtures of the type discussed is not known. 
Some recent papers (Niesel, 1952, Smith, 1956, Reynolds and Hough, 1957) concern- 
ing the mixture formulas have not contributed essentially to the practical solution of 
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the problem here considered. Fig. 10 shows the calculated and experimental permitti- 
vity-temperature dependences of the “BM” system. It seems that the flattening, of 
the experimental responses is due to the fotmation of intermediate boundary layers. 
These may arise in different quantity and composition depending’ on the diffusion 
rate, i. e. the sintering régime employed. 
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Fig. 10. € versus temperature for heterogeneous system “BG“ as measured and calculated 


The measurements carried out show that it is possible to obtain convenient per- 
mittivity-temperature responses over various working temperature ranges by modify- 
ing the composition and sintering régime of such heterogeneous systems, with rather 
high values of the dielectric permittivity. 

Applying the sintering conditions described above, material of a-rather small 
average grain size of about 1—5 microns is obtained. The materials investigated reveal 
closely similar behaviour to that of simple solid solutions of the BaTiO, — SrTiO; 
type. Thus, the samples are ferroelectric with a spontaneous polarization of the order 
of 2 uC/cm?, the coercive force at 50 Hz amounting to about 2 kV/cm. 

The dielectric permittivity versus electrical field responses were obtained from 
HT bridge measurements, at 50 Hz, and showed a typical increase of the permittivity 
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with the field strength, attaining saturation at about 2 kV/cm, when the process of 
domain orientation is accomplished. The biasing field influence, measured by a 1 MHz 
resonance method, showed a typical delay-effect as described for simple solid solu- 
tions. These results seem reasonable, since the material tested consisted of grains, 
each of which is a wimple solid solution revealing the known delay-effect. 


Homogeneous systems 


The material investigated presented a three component system composed of 
barium metatitanate BaTiO,, magnesium metatitanate MgTiO, and barium stannate 
BaSnO,, to be denoted for simplicity as BaTiO, — MgSnO . This system seems not 
to have been investigated in detail, and was chosen since the author expected that 
such a composition might satisfy the requirements as to high dielectric permittivity 
and a low temperature coefficient. Moreover, such materials should not present 
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Fig. 11. Composition of BaTiO,-MgSnO, systems investigated 


a very high power factor. Furthemore, it seemed very interesting to investigate on the 
ferroelectric nature of such systems, a question which should be decided by dielectric 
hysteresis. 

Fig. 11 shows the choice of concentrations pursued in studying the various com- 
positions of the BaTiO, — MgSnO, system, from zero up to 30 mol per cent of MgSnO, 
in BaTiO,. Higher concentrations were not investigated owing to the difficulties of 
producing such samples, which required higher sintering temperatures than could 


Dielectric Investigation of Perosvkite Type Ferroelectrics 405 


be obtained in the actual conditions. Samples of concentrations between 6 and 7 per 
cent, not all of which are shown on the graph, were also produced and measured in 
order to clarify the shape of the Curie-Weiss temperature-concentration curve. 

It is known that barium metatitanate and barium metastannate form ferroelectric 
solid solutions similar to the BaTiO, — SrTiO, system (Smolenski et al., 1951). Simila 
rity in the sizes of the tetravalent Ti and Sn ions (the Goldschmidt ionic radii for the 
coordination number 6 being 0.64 À and 0.74 À, respectively) admits of the formation 
of a perovskite structure with Ba ions. On the other hand, barium metatitanate and 
magnesium metatitanate form no solid solutions (the ionic radii of Ba and Mg being 
1.43 A and 0.78 A), but only heterogeneous mixtures (Shelton et al. 1948) with 
decreasing NT permittivity towards higher MgTiO, concentrations. The presence 
of MgTiO, causes no displacement of the transition temperatures of BaTiO,. 

The grain dimensions were estimated at about 1—5 microns on the average for 
most of the solid solutions prepared. The difference in grain size are not considerable 
owing to the same régime of preparation employed for all samples of this system. 
Most of the samples, particularly those of up to 15 mol per cent concentration, exhi- 
bited a glassy and shining breaking and were about 8 degrees hard by the Moss hardness 
scale. The colour of the samples was yellowish. The samples of concentration exce- 
eding 15 mol per cent of MgSnO, in BaTiO, seemed to require somewhat higher 
sintering temperatures to obtain a low porosity and, consequently, a higher effective 
value of the dielectric permittivity. The breakdown strength for most of the samples 
exceeded 30 kV/cm, which is a rather high value as compared to those given in Table 
I for BaTiO, — SrTiO, solid solutions. 

The dielectric permittivity as a function of the temperature for BaTiO, — MgSnO; 
systems of varying composition, together with the TC values computed from 
these measurements, is shown in Fig. 12. The weak measuring field was of the order 
of 10 V/em, at a frequency of 1 MHz. The samples were maintained in the special 
vacuum holder mentioned above, making it possible to carry out the measurements 
within the wide range of temperatures from —180° up to +180°C. The temperature 
of the sample was measured using a copper-constantan thermocouple whose one 
junction was kept at 0?C and the other was fixed on the face of the sample. The EMF 
of the thermocouple was measured by a compensation method using a Diesselhorst 
compensator, at an accuracy of the temperature measurement of at least 0. EAR 
responses show clearly that what we are dealing with is not a heterogeneous mixture, 
but a solid solution, since the peak value of the dielectric permittivity displaces towards 
lower temperatures with rising MgSnO, concentrations. At the same time the TC e 
values are small for these materials, as the e (T) responses are flattened down at the 
Curie point. 

Thus, from these investigations it is possible to obtain the phase diagram for 
the material tested. As dielectric hysteresis measurements show, the temperature 
corresponding to the peak value of the dielectric permittivity was found to be that 
of the transition from the paraelectric to the ferroelectric phase and thus represents 
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the Curie point of the sample. However, investigations on the ageing effect (Pajak 
1958) show that the Curie point does not represent an intrinsic constant of the 
but is dependent on the actual age of the sample. Thus, instead of the Curie point T., 
as defined by the maximum of e, it seems reasonable to plot the Curie-Weiss Tape 
ture 0 defined as that of the Mosotti catastrophe, and independent of ageing. Fig. 13 
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Fig. 13. Curie-Weiss temperature versus MgSnO, concentration in BaTiO, 


shows the Curie-Weiss temperature versus the concentration for the samples measured. 


One can see that this dependence is nearly linear for concentrations exceeding 2 mol 
per cent of MgSnO, in BaTiO,. For smaller concentrations, the Curie-Weiss tempera- 
ture falls more rapidly with the foreign ions contents. May be this fact is related 
to the changes occurring in the structure of the BaTiO; lattice. The X-rays investi- 
gations now carried out in this Laboratory should decide the problem. 

The effect of biasing on the samples was investigated at room temperature. The 
relative change of the diélectric permittivity, À ele, is not only a function of the biasing 
field and temperature, but is found to be a time-dependent quantity as in the case 
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of previously investigated ferroelectric systems, such as BaTiO, — SrTiO, solid 
solutions (Piekara and Pajak, 1953). In Fig. 14 are plotted the values of A ele versus 
the time, for two solid solutions of different concentrations, after applying and remo- 
ving a biasing field of 10 kV/cm at room temperature. One can see that on applying 
the field the instantaneous increase A & of the permittivity is followed by a slower 
decrease, and that the same occurs on removing the field. The rapid changes of the 
permittivity occur within a few minutes. The character of this effect changes with 
the temperature, and at T>® the quantity A e changes its sign (if the short time 
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Fig. 14. Delay-effects in BaTiO; — MgSnO, solid solutions below and above the Curie point 
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effects in the milliseconds range be neglected) on applying the biasing field. It seems 
reasonable to admit that this behaviour in ferroelectrics of the type considered is 
connected chiefly with an orientation of the domain structure existing below the Curie 
point and, as reported earlier (Anliker et al., 1954), maybe also to some extent above 
the transition temperature. Moreover, the ageing effect described in Part II seems 
to indicate that the delay effects should be related to a domain structure, since 
they are different for young and aged samples. 

Measurements of the dielectric permittivity versus 50 Hz AC field strength for- 
various solid solutions at room temperature were carried out using the HT Schering 
bridge (cf. Part II, Fig. 4.). 

On approaching the Curie point of the solid solution, the following change in 
the shape of the e (E) curve is observed: the maximum corresponding to the beginning 
of the saturation of domain orientation shifts towards lower values of the field strength. 
The rising part of the curves corresponds to the domain orientation; the latter involves _ 
a rise in the polarization, which in turn raises the dielectric permittivity. When satu- 
ration of domain orientation has been attained, any further increase in the field strength 
must appear as a decrease of the permittivity, since the polarisation does not change 
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appreciably. In a true paraelectric, the permittivity should be independent of the 
field strength. However, the materials investigated are not linear dielectrics above 
the transition temperature, and thus a fall in the e (E) curve above the Curie 
point is observed, involved by nonlinearity of the dipole polarization of the ferroele- 
ctric. Theoretical considerations on this problem, for the simple BaTiO, structure, 
are due to Piekara (1954). 

The shape of the e (E) curves was also investigated for a 5 mol per cent solid 
solution at different temperatures (together with Stankowski and Szymanska). It is 
not surprising that this set of curves closely resembles the one described above wherein 
the concentration plays the róle of parameter, at constant temperature. This is explic- 
able since the spontaneous polarization-temperature dependence is similar throughout 
the range of concentrations. 

The BaTiO; — MgSnO, system within the foregoing concentration range proved 
to be ferroelectric, since hysteresis loops appeared for all samples investigated below 
their Curie points. Above this transition temperature, nonlinearity of the polarization 
was found within a large temperature range. Hysteresis loops were observed for 
samples of various concentration at a field strength of 10 kV/cm and a frequency 
of 50 Hz. The oscillograms obtained at room temperature are shown in Fig. 3, Part II. 
The character of the loops changes gradually with the concentration, similarly as 
when varying the temperature of a sample of given concentration. For samples with 
Curie points below room temperature, i. e. of more than 7 mol per cent MgSnO,, 
one still observes some spontaneous polarization; at higher concentrations nonlinearity 
of the dielectric polarization is present, vanishing at concentrations exceeding 20 mol 
per cent. The room temperature values of the total P and spontaneous polarization P, 
as well as those of the coercive force E,, which also decrease with the MgSnO, contri- 
bution, are collected in Table II. It is interesting to note that the values of the spon- 
taneous polarization (estimated graphically by tracing the tangent at the point cor- 
responding to saturation of the domain orientation (Stankowski, 1958)) are relatively 
high. This means that the. system possesses a markedly ferroelectric character similarly 
to the other perovskite type ferroelectrics. 

The temperature dependence of the total polarization P and the spontaneous 
polarization P, for a quite young (see Part II) 5 mol per cent sample is shown in Fig. 15. 
It is clearly seen that the spontaneous polatization P, vanishes some centigrades 
above the temperature of the maximum value of e, which is 50°C. The total polari- 
zation P diminishes more rapidly in the paraelectric than in the ferroelectric phase. 

The TC e values are sufficiently small to present the possibility of applying 
BaTiO, — MgSnO, solid solutions for technical purposes e. g. as condenser materials. 
A comparison of the BaTiO, — MgSnO, solid solutions and the BaTiO, — SrTiO; 
systems proves that the former constitute a class of ceramics presenting an advan- 
tageous TC e, since between — 180°C and +180°C their |TC e| is smaller than 
200 x10-4 deg-1. The dielectric permittivity of the MgSnO, — BaTiO, systems 1s 
rather high, exceeding 4000 for young samples and 2000 for aged ones. 
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From a theoretical point of view it is difficult to explain the low TC e value 
near the transition temperature as obtained by replacing the Ti and Ba ions of the 
BaTiO, crystal structure with Sn and Mg ions. respectively. As already mentioned, 


there is as yet no exact theory able to account for all the ferroelectric properties even 


TCe, = TCez = TCe, 


SI 
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_ Fig. 16a 


of pure BaTiO,. The present problem is more complicated because of the lattice 
distortions introduced by the Mg and Sn ions. For a ferroelectric obeying the Curie- 
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Weiss law € = og the paraelectric phase, the definition TA een leads 
to the hyberbolic TC e (T) dependence 
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This means that the TC e values are independent of the C value characterizing the 
ferroelectric (Fig. 16a). The Curie constant C is decisive for the e value and it is clear 
that the maximum e/TC e ratio can be obtained only if the value of the constant (H 
is large. Similarly, for a ferroelectric obeying the Curie-Weiss law i 
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Fig. 16b 


Fig. 16. € and TC e versus temperature calculated for the paraelectric and ferroelectric phases 


in the ferroelectric phase, the relation TC e = 


5 T holds. Since the O's for both 
phases are different, the TC e (T) curve is obtained as in Fig. 16 b. Thus, beyond 
the Curie-Weiss region, the TC & values are the same for all ferroelectrics for the 
same @—T values, and advantegous e/TC e values are determined by the constant C 
only. Within the Curie-Weiss region, at the transition temperature T,, the TGe 
changes its sign suddenly (dashed line in Fig. 16 b). The maximum TC e values are 
thus determined by T,. Such experimental TC e(T) curves are observed in BaTiO, 
and solutions of the BaTiO,— SrTiO, group, and are shown in Figs. 4 and 12. In 
ferroelectric systems such as BaTiO, — MgSnO,, for which the difference 0, — 0, is 
considerable, the values of TCe with 0, — T < 0, are small and thus the TC e (T) — 
curve is flattened. 

A possible suggestion as to another cause of the small TC e values observed 
in BaTiO, — MgSnO, systems is that it resides with the character of the transition, 
which must be decisive for the TC e values in the neighbourhood of the transition | 
temperature. The transition is situated within a range of sharp e (T) dependence 
involving high TC e values. Since the ceramic material consists of grains which may 
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somewhat differ in composition, therefore for every particular grain the transition 
temperature may also be somewhat different. Thus, the sample as a whole may present 
a transition extending over a range of temperatures, leading to greater deviation from 
the Curie-Weiss law, involving the flattening of the e (T) dependence in the vicinity of 
the actual Curie point of the sample. Thus, the character of the transition is now 
somewhat different and the TC e values observed are lowered. From the point of 
view of crystal structure, such an interpretation seems reasonable, since the foreign 
ions — chiefly the Mg ions, which are considerably smaller that the Ba ions — do not 
favour coupling of the ferroelectric dipoles, which is the cause of the spontaneous 
polarization below the Curie point. Thus, the latter increases with decreasing tempe- 
rature less rapidly than in the case of pure BaTiO,. This is indicative of the fact that 
we are dealing with a type intermediate between continous and discontinous transition. 
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DIELECTRIC INVESTIGATION OF PEROVSKITE TYPE 
FERROELECTRICS 


Part II: Ageing Process in Ferroelectrics 
By Zpzistaw PAJĄK 


Polish Academy of Sciences, Institute of Physics, Dielectric Laboratory, Poznan 
(Recived February 16, 1959) 


The dielectric behaviour of BaTiO, type ferroelectrics in the process of ageing 
was investigated and new ageing effects were detected. It was found that the Curie point 
(i. e. the temperature of the paraelectric-ferroelectric transition) is displaced towards 
higher temperatures and the spontaneous polarization, coercive force and hysteresis 
losses vanish in the process of ageing. 

A domain mechanism of ageing leading to the formation of domain antiferroelectrics 
is proposed, and some related effects such as rejuvenation and delay effects are discussed. 


Introduction 


Ageing is the term used to denote the process which develops spontaneously 
with time following the preparation of a dielectric and which leads to the state of 
lowest free energy. In the course of ageing, time variations of certain physical pro- 
perties of the dielectric, such as the dielectric polarization, dielectric permittivity, 
power factor and others, are observed. Notwithstanding the fact that ageing effects 
are of great importance to the practical application of ferroelectric materials, few 
experimental and theoretical publications on the problem are in existence. It is known 
that the dielectric permittivity and power factor decrease, subsequent to the prepara- 
tion of BaTiO, ceramic samples (Kazarnovski, 1952, Piekara and Pajak, 1952, Novo- 
siltsev, Khodakov and Shulman, 1952, Bogoroditski and Verbitskaya, 1953). From 
recent investigations (Plessner, 1956), this appears to be a logarithmic process. 
Moreover, ageing of such properties as the piezoelectric constant da and the electro- 
mechanical coupling factor was reported for some ferroelectric materials (Mason, 1955). 


Curie point and Curie-Weiss temperature in ageing process ` 


It was found by the present author that the dielectric permittivity and power 
factor of BaTiO, — MgSnO, solid solutions also undergo ageing. To observe the 
changes over a greater range of temperatures, measurements of the dielectric permit- 
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tivity as a function of the temperature were carried out for young samples and ones 
aged at room temperature. Measurements were made when heating the sample to 
avoid temperature hysteresis of the dielectric permittivity (Piekara and Pajak, 1952). 
The results of such measurements are shown in Fig. 1. It is seen that the ageing effect 
is limited to the ferroelectric phase of the material tested. From Fig. 1, a new effect 
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Fig. 1. e versus temperature for a young and aged BaTiO, — MgSnO, solid solution S-O4 (for concen- 
tration, see Table II, Part I) H 


reported earlier (Pajak, 1958) is seen to exist, consisting in a displacement A T, of 
the Curie point towards higher temperatures due to the ageing of the sample. By 
the curves, this experimental shift in the Curie point is considerable and amounts 
to about 20? centigrade. Thus, the Curie point i. e. the temperature of the ferroelec- 
tric-paraelectric transition is also a property subject to ageing. It should be stressed 
that the Curie point, similarly to the dielectric permittivity, power factor and mecha- 
nical coupling factor, does not constitute an intrinsic property of the material, but 
depends on the history of the sample and is cooperative in describing its actual state. 
Thus, the common use of the Curie point as a quantity representing the temperature 
of thermodynamical equilibrium between the paraelectric and ferroelectric phases 
in the phase diagram is disadvantageous in that, for different times following produc- 
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tion of the solid solutions, the diagrams have different shapes. On the other hand, 
on comparing the paraelectric part of the e (T) curve in Fig. 1,.one can see that ageing 
does not occur in this phase; this means that the effect is related to the ferroelectric 
properties only. It follows that the Curie-Weiss temperature 6, estimated by extra- 
polation of the 1/e (T) curve in Fig. 1 remains constant during the process of ageing. 
Thus, it seems reasonable to admit this latter quantity as a true constant characterizing 
the material better than the Curie point T.. 

It seems possible to interprete the process of displacement of the Curie point 
in the course of ageing as follows. The Curie point determined experimentally as 
the temperature of the peak value of the dielectric permittivity is the one corresponding 
to the intersection of the Le (T) curves of the ferroelectric and paraelectric phases 
(Fig. 2). Since the 1/e (T) curve of the paraelectric phase, and therefore the Curie- 
Weiss temperature 0,, do not vary in the process of ageing, any change in the Curie 
point T, must necessarily result from a change in the Le (T) curve of the ferroelectric 
phase. It follows from. the measurements that, owing to the lowering of dielectric 
permittivity in the aged ferroelectric phase, the 1/e (T) curve generally shifts parallel 
towards higher temperatures. As a result, the temperature of intersection of the 
two 1/e (T) curves (T°) shifts towards higher temperatures (OT). What, then, is the 
cause of the change in the 1/e (T) curve for the ferroelectric phase? It seems that 
it may reside in a spontaneous lowering of the P, value for the sample as a whole 
due to domain compensation, the elementary cell retaining its dipole moment un- 
changed. The reason for which such an assumption can be made is that a decrease 
in P, during ageing has been found experimentally (see below). Such a lowering of 
the P, value should correspond to a shift of the Le (T) curve toward higher tempera- 
tures, as when removing hydrostatic pressure (Klimowski and Pietrzak). Thus, the 
ageing process seems to consist in a spontaneous release of the mechanical stresses 
stored chiefly in the process of passing through the transition point to lower tempera- 
tures, when the spontaneous polarization appears. Therefore the Curie point T, for 
young samples is lowered owing to the stresses involved, similarly as in the case of 
hydrostatic pressure acting in the sample (Merz, 1950). Ageing corresponds to the 
removal of pressure involving displacement of the Curie point towards higher tempe- 
ratures. | 

The fact that the transition appears at different temperatures seems to be expli- 
cable also from the point of view of the theory of spontaneous polarization in barium 
metatitanate proposed by Piekara (1958). Although this theory is given for a one- 
domain BaTiO, single crystal, the general conclusions hold for other perovskite type 
_ ferroelectrics. To give an exact theoretical explanation of the displacement of T,, the 
theory of polarisation should be generalized for multidomain crystals ov, still better, . 
for polycrystalline material. However, no such theory exists at present and, for a first 
qualitative explanation, the theorical results for one-domain crystals will be used. 

Within the Curie-Weiss region (between both catastrophe temperatures 0, and 
0, of the paraelectric and ferroelectric phases) the P, (T) function is not monovalent, 
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SO ————— 


which means that P, can appear or vanish at any temperature within this region 


(Fig. 2), i. e. two states, namely P, > 0 or P, = 0, are possible. The probability of 
the transition from the state with P, > 0 to that with P, — O at a given temperature 
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produce a rather large value of E,, whilst in an aged sample Æ} vanishes. Thus, in 
a young sample, heating to the Curie-Weiss region (0, 0;) produces the transition 
P, 0 P, = 0 at a lower temperature than in the case of an aged sample. The 
extremal value for the ferroelectric transition temperature in the quite young sample 
seems to be 0,. Experiments performed with ageing BaTiO, ceramic proved that 6, 
(determined aoe the l/e (T) dependence) remains constant, while 9, shifts SSC 
towards higher temperatures (0,— Oj). At the same time, the transition tempera- 
ture T, (determined as the temperature of the peak value of e) shifts towards 6” 


! 
Ce Le) 
Dielectric hysteresis in ageing process 


Investigations on the dielectric hysteresis loops for young and aged samples 
were undertaken to clarify the behaviour of the spontaneous polarization in the pro- 
cess of ageing. The results of observations made together with J. Stankowski have 
been reported recently (1958). Fig. 3 shows a set of oscillograms obtained with various 
samples at room temperature previous and subsequent to ageing during one year. 
It is clearly seen that ageing results in the vanishing of the spontaneous polarization 
PY. At the same time one can observe the vanishing of the coercive force E, and hyste- 
resis losses 4. As ageing proceeds, the hysteresis loop becoms constricted, gradually 
assumes the shape of a double loop and, later on, that of a one-branch curve, still 
showing nonlinearity of the polarization; finally, the P (E) dependence becomes 


- linear. This behaviour was observed also for simple BaTiO,, for BaTiO,-SrTiO, solid 


solutions and heterogeneous systems and probably is characteristic of other perov- 
skite type ferroelectrics. 

The dielectric hysteresis investigations were completed by HT bridge measure- 
ments of the dielectric permittivity as a function of strong electric fields for young 
and aged samples. Fig. 4 brings two examples of e (E) curves for pure BaTiO, and 
for a 5 mol per cent solid solution of MgSnO, in BaTiO . One can see that for a young 
sample the dependence is strong. In an aged one, the maximum of the curve shifts 
towards higher field strengths, which means that saturation of the domain orientation 


requires more energy than in the case of a young sample. These results are in good 


agreement with the fact that the movement of a 180? wall separating compensated 
domains requires more energy than the movement of a 90° wall separating non- 


-compensated domains. 
Mechanism of ageing 


The results of the present investigations serve to confirm the following conception 
of the mechanism of ageing. Primarily, the spontaneous polarization P, falls, or, in the 


- extreme case, vanishes; similarly, the coercive force E, and hysteresis losses A vanish. 


NEED ur E s 
1) Similar ideas were put forward by Bogoroditski and Verbitskaya (1953), with no experimental 


proof however. 
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Fig. 3c. Hysteresis loop for young and aged BaTiO, — MgSnO, solid solutions (for concentrations, 
see Table IT, Part D 


Dielectric Investigation of Perovskite Type Ferroelectrics 515 


4000 


2000 


EkV/km 


| 0 Í 2 -à 4——EkV/cm 
Fig. 4b. & versus 50 Hz measuring field strength for young and aged BaTiO, and S-05 solid solution 


As explained theoretically by Mason (1955), the vanishing of P, would necessarily invo- 
lve a decrease in the dielectric permittivity, power factor, and electromechanical coup- 
ling factor, these being secondary effects of the ageing process. The state with a high 
value of P. as obtained after passing through the transition point to lower tempera- 
tures, may be said to represent a non-stable state possessing additional energy with 
respect to the state P, — 0. At the same time, the mechanical strains are stored in 
the sample owing to tetragonalisation of the initial cubic structure, as well as to the 
presence of the chiefly 90° domain walls possessing wall energy. Moreover, the depola- 
rizing field E, resulting in the sample assumes a large value and thus favours the 
change of the P, direction of the domain. It is known from optical observations that 
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the domain walls shift spontaneously when the state of energy existing in the crystal 
is disturbed in dny way. The measurements described above show that spontaneous 
wall motion takes such a direction as to reduce the P, value and mechanical strains 
by making the 90° domain walls vanish. It is possible to realise an antiferroelectric 
or otherwise compensated domain arrangement so as to obtain P, — 0 for the sample 
as a whole. A crystal thus aged may be termed a domain antiferroelectric. Simultane- 
ously, the E; value decreases as a result of the decrease in the number of uncom- 
pensated domains. Hence, the hysteresis losses diminish because the majority of the 
domains are excluded from the process of reorientation under the action of the applied 


field. 


Rejuvenation 


. Ageing appears to be a reversible process, since the stable state of domain arran- 
gement may be disturbed by certain external factors acting on the aged sample. Such 
disordering factors are, in the first place, the temperature and external electric field. 
[t may be that mechanical stresses also play some part in the process. The latter may 
be conveniently termed rejuvenation, as consisting in a reversal of the ageing process. 
Rejuvenation leads to a reincrease in P,, E, and A, which, in turn, involves-analogous 
changes in the other properties, such as the dielectric permittivity, power factor 
and so on. Thus, the experiments performed show that the low energy state (aged 
sample) can be raised by introducing thermal or electrical energy, which dissipates 
to destroy the coupling of the domain ensembles formed in ageing. To this purpose, 
thermal rejuvenation seems to be the most convenient, since it acts isotropically, 
whereas the AC field acts anisotropically favouring an unidirectional arrangement 
of the dipole moments throughout the sample. Thus, rejuvenation under the influence 
of the field is not so complete as that brought about by high temperatures. The 
temperature and time required to bring about rejuvenation depend on the concen- 
tration of the solid solution; however, in most cases, heating to about 100°C above 
the Curie point for ca 15 minutes is sufficient. 

From the standpoint of these experiments, the thermal hysteresis of the dielectric 
permittivity as observed earlier by Piekara and Pajak (1952), now seems to present 
a particular case of the rejuvenation process. Heating the sample above the actual 
Curie point causes a reversal of the ageing process, resulting in a replacement of the 
dielectric permittivity and power factor and in a replacement of the Curie point 
towards lower temperatures i. e. in the direction opposite to that of the displacement 
in the ageing process. The experimental results, shown in Fig. 5, confirm this state- 
ment. The thermal hysteresis of the permittivity was observed for one and the same 
sample of 6 mol per cent MgSnO, in BaTiO, previous and subsequent to ageing. 
One can see that A T, after one year of ageing is approximatoly 14°C, whereas for 
a young sample its value amounts to only about 2°C, just as the difference between 
the e (T) curves from the heating and cooling of an aged sample is considerably grea- 
ter than that obtained with a young sample. j 
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Fig. 5. Thermal hysteresis for a young and aged S-06 solid solution 
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The Curie point replacement towards lower temperatures observed in the thermal 
hysteresis of € results from restoring of mechanical stresses involved by the process 
of rejuvenation. 

Moreover, it is interesting to note that the ageing process is more pronounced 
in the distorted BaTiO, structures (solid solutions), whereas rejuvenation is easier 
in pure BaTiO. This means that domain compensation or coupling is more marked 
in the presence of lattice imperfections, E. g., the behaviour of a BaTiO, sample 
and of a5 mol per cent MgSnO, solid solution in BaTiO, will be compared: a 10 kV/cm 
AC field of 50 Hz applied during 10 minutes rejuvenates the BaTiO, sample quite 
distinctly, restituting the normal hysteresis loop, whilst in the solid solution reju- 
venation leads no further than to the reappearance of a double hysteresis loop. In 
the case of BaTiO, an interpretation of the rejuvenation process may be suggested 
by which the external energy is used for the activation of the domain. The replace- 
ment of the original ions by Mg and Sn must necessarily modify the internal potential 
of the lattice. If this leads to an increase in the potential barriers of the domain walls 
in the aged sample, rejuvenation in the case of solid solutions requires more energy. 


Delay effects as an ageing process 


From the point of view of such a mechanism of ageing, it seems possible to suggest 
the following new explanation of the previously observed delay effects i. e. the time 
changes of the dielectric permittivity under the influence of biasing fields (Piekara 
and Pajak, 1953). These may now be treated as a short time ageing effect. Measure- 
ments of the relative change of the dielectric permittivity A ele as a function of the 
time, on applying or removing the biasing field show that these responses depend 
on the state of ageing of the sample. Fig. 6 shows the responses at room temperature 
for BaTiO, and for one of the heterogeneous mixtures. These measurements were 
carried out at 1 MHz frequency, the biasing field being 5 kV/cm, the measuring 
field — 1 V/cm. The bias was applied during 5 minutes, and after its removal the 

ermittivity was measured during the same period of time. Measurements were 
carried out for both directions of the field; the differences, however, are not essential 
ones, and for clarity the results are plotted only for one direction. One can see that 
at strong fields the & values are smaller for all samples when aged. This fact seems 
comprehensible if it be assumed that the delay effect is one of disturbing the arrange- 
ment of the domains by the biasing field involving a change in the P, value. Appli- 
cation of the field causes a rearrangement of the domains whose P, directions differ 
from those of the field. This bias can temporarily create a number of domains whose P, 
direction is perpendicular to the E direction, which involves an increase in the permit- 
tivity value, as it is known that e is much greater in the c-axis direction than in 
that-of the a-axis. Such an effect takes place in the first phase of the field action 
or when the field strength is not too high. The stresses thus involved by the creation 
of 90? domain walls vanish with time, leading to the second phase of the effect i. e. 
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the slow orientation of the domain whose P, is perpendicular to the field direction, 
whereas e decreases, the e, value being the lowest of the possible average values of e. 
If the field strength is sufficiently high, the first phase proceeds so rapidly (in about 
l millisec, as observed by Kaczmarek) that it cannot be observed by the method 
described, and the second phase i, e. the rearrangement of the domains according 
to the field direction prevails involving a decrease of the & value only. In this case 


the domains cannot be reorientated by the weak measuring field, and hence the de- 
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Fig. 6. Delay effect for young and aged BaTiO, (S-00) and for the “BM” — heterogeneous system 
(see Part I) 


crease in permittivity is generally observed. Thıs decrease proceeds with time as the 
remaining domains switch over in the field direction. In the extreme case of high 
biasing voltages, saturation of the permittivity decrease occurs, i. e. the minimum 
value of g is observed a very short time after the bias is applied and this value does 
not vary further, as has been confirmed by numerous experiments. 

The effects observed on removing the field can be explained by the same mecha- 
nism of domain orientation accompanied by a vanishing of the stresses. On removal 
of the bias, the domain returns to the low energy arrangement (antiparallel or other- 
wise compensated state). Besides immediate switching, this reorientation may also 
be attained through a P, oriented perpendicularly to the measuring field (90° domain) 
involving a transient increase of the permittivity, followed by the subsequent decrease 
resulting from ageing. The results shown in Fig. 6 are to be explained by accounting 
for the fact that, in the aged sample, one is dealing with compensated domains, and 
thüs the field applied is unable to switch over so many domains as in a young sample. 
Namely, in the aged sample the antiparallel coupling (180? domains) must first be des- 
troyed and the domain may go over to the field direction by way of the perpendicular 
orientation. Hence, for equal times following application of the bias, it is the aged 
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sample that exhibits the higher e value owing to the existence of 90° domains. It was 
found that, for sufficiently strong fields, these differences vanish because most of 
the domains become aligned to the field direction causing the decrease of e. Thus, 
the delay effects may be considered to represent a particular case of short time ageing 
effect. 
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SUR LA TENSION SUPERFICIELLE DE L’ALCOOL ETHYLIQUE 


M. Bonwras ET I. BABUTIA 


Laboratoire de Physique; Institut Polytechnique, Timisoara 


(Recu le 8 mai 1959) 


On montre que la tension superficielle de l'alcool éthylique mis en rotation est 
plus grande qu'en repos et qu'elle n'est pas une fonction linéaire de température. 


On a montré dans quelques travaux antérieurs b?9* que la tension super- 
ficielle des liquides change s'ils sont mis en rotation. Elle a, pour la plupart des cas, 
des valeurs cinétiques plus grandes que les valeurs statiques. Cependant, dans cer- 
tains domaines de température, concernant quelques substances; comme par exemple 
l'eau, ou le glycol, quand le liquide se trouve en rotation, les valeurs de la tension 
superficielle deviennent moindres que les valeurs en repos?. 

Dans ce travail nous avons étudié la tension superficielle de l'alcool éthylique 
en rotation, dans le domaine de température de 6,5°C à 13°C. Parce que ce sont les 
valeurs en rotation comparées aux valeurs statiques qui nous intéressent, nous n'avons 
déterminé que les angles à la rupture, par la méthode Lecomte du Noüy, utilisée 
dans les travaux cités. L’appareil a été adapté en vue de mesurer la température du - 
liquide tout prés de l’anneau*. De méme, on a eu soin d'assurer l'homogénéité élas- 
tique du fil de torsion?. 

La figure représente les valeurs de l'angle à la rupture, mesurées alternativement 
en repos (e) et en rotation (O). + | 

On voit que dans le cas des mesures où le liquide se trouve en repos, la dépen- 
dance entre l’angle et la température est rectiligne, correspondant ainsi aux données 
de la literature®. Au contraire, dans le cas des mesures effectuées durant la rotation, 
la courbe des valeurs présente un minimum autour de la température de 10°C, en 
s'approchant de la droite des valeurs statiques. 

Il est évident que la formule de Eotvos-Ramsay-Shields ne peut pas s'appliquer 
quand le liquide se trouve en rotation. Cela parait ne pas étre une anomalie, mais 
une propriété d'un grand nombre de liquides pures, sinon de tous. Des recherches 
futures auront comme objet d'élaircir ce probléme. 

Nous remercions le professeur A. Cisman pour l'aide matériel de grande valeur 
qu'il nous a accordé. 
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VISCOSITY OF LIQUID HYDROGEN SULPHIDE 


By J. W. HENNEL AND K. Krynicxki 
Cracow Centre of Nuclear Physics, Institute of Nuclear Research 
(Received May 2, 1959) 


The viscosity of liquid H,S under the pressure of its saturated vapour in the tempe- 
rature region from — 11.5 to 50°C was measured. A closed type, pressure resistant Ostwald 
viscometer was used. 


Hydrogen sulphide is a very well examined chemical compound; nevertheless, 
the viscosity of its liquid phase has not been known for temperatures above — 63°C. 
The only measurements were carried out by Steele et al. (1906) in the temperature 
region between — 80° and — 63°C. 

Hydrogen sulphide is especially suitable for investigations of nuclear spin relaxa- 
tion phenomena. This is due to the simplicity of the theory in the case of molecules 
containing two protons. According to the theory of Bloembergen, Purcell and Pound 
(1948), relaxation processes in liquids are close related to viscosity. This makes 
the viscosity of liquid H,S a matter of considerable interest, especially in the region 
near room temperature where paramagnetic nuclear resonance experiments are easier. 


, Evolution of Hydrogen Sulphide 


Hydrogen sulphide was obtained by the addition of a 15 per cent HCI solution 
to an aqueous solution of chemically pure Na,S. The latter solution was freed from 
SO,+* and COg4** ions by treatment with excess BaCl, and subsequent filtering. 
Evolving hydrogen sulphide was made to pass through three washers containing 
water, Na,S solution to remove HCl, and water, respectively. 

It was then dried in a spiral glass tube immersed in a cooling mixture consisting 
of solid CO, with alcohol. Dry H,S was collected in a sealed glass container cooled 
by liquid air. The air above the solid H,S was removed with a rotary oil vacuum 
pump. H,S was subsequently distilled into the viscometer. 


Description of the Method 


A closed type Ostwald viscometer was used. This viscometer. which is shown 
in Fig. 1, resembled the one described by Lewis (1925). To ensure that the viscometer 
would withstand the pressure of saturated H,S vapour at 50°C (35.6 Atm.), no glass 
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bulbs were used and all parts were made from tubings. Tubes a and b were of 7.3 mm 


outer and 5.3 mm inner diameter, and tube c — of 5.7 mm outer diameter; d is the 


capillary tube 8 cm long and of 0.39 mm inner diameter. The total length from c 


a.m 


CIN: 


[^ 
Fig. 1 


length of the 


to b was 28 cm. On the viscometer, two marks were made by means 
of nitrocellulose paint (u.m. — upper mark, l.m. — lower mark). 

While the hydrogen sulphide was vacuum distilled from the main 
container to the viscometer, the former remained at a temperature 
of about — 20°C, and the viscometer was immersed in a mixture of 
solid CO, and ethyl alcohol. After the viscometer had been filled 
with the required amount of liquid, both the container and the 
viscometer were quickly cooled by liquid air. Then the viscometer 
was sealed off at point e, yielding a closed system containing air 
free HS, 

At room temperature the viscometer was inserted into a thick 
walled glass tube of 23 mm inner diameter and 42 cm long, through 
which water from the thermostat was allowed to flow. This tube was 
mounted in a suitale bohlder. For the security of the observer, 
a perspex plate 1 cm thick was inserted between him and the whole 
apparatus. 

To measure the viscosity, the viscometer with the outer tube was 
inverted, allowing the liquid to rise in tube b until it was nearly 
completely filled. The viscometer was then replaced right, care being 
taken that the space c-d was completely filled with liquid, and the 
time required for the liquid to flow through the tube between the 
upper mark and capillary was measured. This time ranged from 150 sec. 
to 280 sec. for water and from 25 sec. to 32 sec. for H,S. A calibrated 
stopclock was used. The temperature was measured with an accuracy 
of + 0.1°C by means of a calibrated glass-mercury thermometer 
placed within the outer tube. 

The formula used for computing the viscosity was the one propo- 
sed by Ostwald. When the kinetic energy corretions are made, the 
formula cau be written as follows (Lewis 1925, Kohlrausch 1955): 


n = (0, — 0,) tk — mo, V/8x tl 


wherein g,, o, are the densities of the liquid and the vapour, res- 
pectively, ¢ is the time of flow, V is the volume of flow, J is the 
cappillary in cm, k is the viscometer constant, and n is the visco- 


sity of the liquid. 
For the constant m, a value of 0.5 was assumed because of the cone-shaped ends 
. of the cappilary tube (Kohlrausch 1955, p. 158). To verify this value and to estimate 


the systematic 


errors from other sources, the viscometer was filled. with carefully 


purified ethyl ether by the method described above. The results obtained for the 


hy ee - 
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viscosity of ethyl ether shown in Table I are in good agreement with values taken 


from Landolt-Bórnstein (1955) p. 198. The differences amount to about 2 per cent 
and less. 


TABLE I 
Control measurements. Viscosity of ethyl ether (cp). 


Authors’ 


°C A Landolt 
measurements 

10 0.2582 0.2602 

20 0.2324 0.2360 

40 0.1935 0.1970 

70 0.1496 0,1528 


The constant k depends on the volume of liquid contained by the viscometer. 
As the volume of the liquid phase in the sealed viscometer changes with the tempe- 
rature, the viscometer was calibrated using different volumes of carefully distilled 
water. This was done at a temperature of 20°C, and the viscosity of water was assumed 
to be 1.005 cp (Hodgman 1956 p. 2030). The volume of the liquid actually contained 
by the viscometer was expressed in arbitrary units by measuring the distance 6 between 
the level of the liquid and the lower mark. As a result of calibration, a plot of E versus 
6 was drawn. : 

The density of liquid H,S was taken from measurements by Baxter et. al (1934). 
The density of saturated H,S vapour was calculated from the Van der Vaals equation. 
The constants a and b were taken from Hodgman (1956, p. 2127). 

The thermal expansion of glass was not taken into account, as its influence on 
measured 7 values does not exceed 0.1 per cent 


Results and Discussion - 


The results obtained by the foregoing method and expressed in centipoise are 
shown in the following table: 


TABLE II 
Viscosity of liquid hydrogen sulphide 


15.0 20.0 30.0 40.0 50.0 °C 
0.133 0.126 0.116 0.105 0.094 | cp 
The estimated error is less than 4 per cent. The main sources of error is the strong 
dependence of the constant k upon the volume of liquid contained, and the rather 


small degree of accuracy in estimating this volume, resulting from the fact that the 
viscometer was made only from tubings without bulbs. 
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The Reynolds number R, = 2Re, w/n, where w is the mean velocity of the 
liquid in the capillary tube and R is the radius of the cappillary, was 2040 and 2650 
for hydrogen sulphide at the temperatures 20°C and 50°C, respectively. In the case 
of ethyl ether, the number was A58 at 20°C. While the liquid was in motion the velocity 
decreased; however, in computing the Reynolds numbers the maximum value of the 
velocity was used. According to Prandtl (1949), the critical value in the case of rounded 
inlet edges is of the order of 40000. 

The authors wish to thank Professor H. Niewodniczariski Ph. D. for rendering 
the present investigation possible at his laboratory, and are indebted to Mr. S. Moroz for 


producing the viscometer. 
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A NEW METHOD OF GAMMA BACKGROUND ERADICATION 


By J. Benisz AND Z. KIERZKOWSKI 
Higher Pedagogical School, Katowice 


(Received April 29, 1958). 


The three most important methods of eradicating the gamma background consist 
in 1) under-development, 2) oxidation of the nuclear emulsion, 3) infrared irradiation. 
In the first of these, the plate is not fully developed. The second is based on the oxidi- 
zing effect of a) water vapour (Albouy, Farage 1951; Violet, Gilbert, Barkas 1951) 
of appropriate temperature acting on the emulsion over a determined period, b) hydro- 
gen peroxide — in this case the plate is placed for some hours over a H,O, solution 
of several 9595. c) chromic acid — the plate is introduced into a 2% solution of the 
acid (Beiser 1952) for an appropriate period. In the last method, Herschel's effect 
is utilized (Jarezyk, Lewandowski 1956). The foregoing methods, however, result 
in a weakening of the tracks of ionizing particles. At this laboratory, a new method 
of eradicating the gamma background using acetic acid has been evolved. Nuclear 
K-2 type plates, 100 u thick, produced by ..Agfa", were impregnated with lithium 
citrate and irradiated with neutrons, yielding triton and & tracks according to the’ 


reaction: l 
on + sli = $He + iT (1) 


Subsequently, the plates were impregnated with a 0.1% acetic acid solution 
for 10, 20, 30, 40, 50, 60, 707 110, 120, 180 minutes and developed with amidol. The 
background was evaluated with a Zeiss grid inserted into the eye-piece of the micros- 
cope, the object lence magnifying 110 times, and the eye-piece 15 times. A conside- 
rable decrease in the gamma background occurred at longer periods of impregnation 
with 0.1% acetic acid. However, the number of grains on the triton and (recoil) proton 
tracks remained unchanged within the limits of error. 

In assessing the results, the signal to background ratio was taken into account 
(Jarczyk, Lewandowski 1956). The signal is defined as the number of grains on the 
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track of an ionizing particle per unit area, the background being the number of grains 
of the background per unit area. As the number of grains on the tracks of ionizing 
particles was independent of the time of impregnation, the ratio of the signal to the 
background may be raplaced by the background alone. Fig. 1 shows a graph of T,/T 
as function of the time of impregnation. 
T, and T denote the number of grains of the background per unit area of a standard 
plate and of one impregnated with 0.1%, acetic acid, respectively. 

It is seen from Fig. 1 that the T,/T ratio increases towards higher values of the 
time of impregnation, attaining saturation at t — 120 minutes (at this point, the 
background is one halt its original value). Acetic acid can be presumed to affect the 


^us 


Fig. 1 


gamma background in a nuclear emulsion either through a) a process of under- 
-development brought about by the low pH value of the impregnating solution, or b) 
through tanning the emulsion. 

To decide which factor is essential in desensitizing the emulsion, some K-2 plates 
were under-developed following irradiation. From the results obtained, the signal 
to background ratio did not exceed 14 at the best. On comparing this result with 
those obtained with acetic acid, it can be concluded that: 


a) the decrease in gamma background is essentially due to a local effect of tanning; 
b) the effect of under-development is a lesser one. 
Experiments carried out by the authors prove that if a nuclear emulsion is impre- 


gnated with an appropriately prepared solution of acetic acid for a sufficiently long 
time, a considerable decrease in the gamma background can be obtained. The number 
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of grains in a proton or triton track is not modified hereby, within experimental 
error. Lie 

The authors wish to express their indebtedness to Mr. L. Jarczyk for his valuable 
remarks, and to Dr K. Grotowski for irradiating the plates with neutrons at the Nuclear 
Research Institute at Cracow. 
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HIGH PERFORMANCE LABORATORY ELECTROMAGNET 
Bv Lupwik Nowicki 
Institute of Physics of Polish Academy of Sciences, Department of Ferromagnetics in Poznan 


( Received April 27, 1959) 


The present note concerns the design of a high performance electromagnet 
suitable for radiospectroscopy work. The magnet consists of a U-type yoke of low- 
carbon steel, Armco poles and water-cooled coils. The device is mounted on a rotable 
support provided with a protractor yielding an accuracy of 0.5°, and with a clamping 
screw. The section of the yoke has an area of 130 cm?, the poles and pole. pieces have 
a diameter of 100 mm, and the gap is adjustable within the range of 10 to 90 mm. 


1 
` 


Fig. 1 Part view of coil body > Ee 


A new and simple system of cooling has been adopted. Flat cooling channels for the 
circulation of water have been cut in the brass body of the coils; the channels, which 
are shown in Fig. 1, ard covered with thin plate. Fig 1 shows a part view of the body — ` 
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of the coils. The solution adopted yields proper cooling of the entire surface notwith- 
standing the restricted wall thickness of the body, amounting to no more than 7—9 mm. 
The construction, which is a simple one, requires less space than the one proposed 
by Piitter (1955), thus making it possible to achieve 2100 windings of 2 mm wire 
in each of the coils. Flat coils and windings of rectangular section have been adopted. 
This involves a somewhat greater resistivity; however, the cooling area is larger. 


Fig. 2 General view of electromagnet 


Moreover, coils of such form present an economy of space and dissipated flux, and 
the gap is more easily accessible. The total dissipated flux amounts to about 7% 
with a 40 mm gap and a field of 6 kOe. In constructions used at present, it is usual 
for the dissipated flux to be 2 to 4 times graeter. The excitting coils act as a solenoid 
magnet giving a homogeneous field of up to 4 kOe in a volume of about 1 1. The 
general view is seen in Fig. 2. The magnetization curves are shown in Fig. 3a; for 
comparison, the recently published magnetization curves of Lange’s electromagnet 
are given in Fig. 3 b. Tables la and 1b bring the data of the electromagnet described 
and for that of Lange, respectively. It can be seen therefrom that the former gives 
stronger fields in considerably greater volumes, notwithstanding the fact that its 
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Fig. 3a Magnetisation curves at different gaps poles diameter @ = 100 mm 


TABLE 1. 


Parameters of electromagnet described in present note (a) as compared to those of Lange electromagnet (b) 


Total number of windings 
Total resistance as of windings in series 


Maximum exciting current 


Maximum supply power 


Maximum excitation 59000 A-windings 44000 A-windings 


Diameter of poles 
Gap 

Total length 
Total height 


Cooling water flow 


Total weight 


a, 4200 b, 2200 
16.7 1l 
14 A 20 A - 
(short duration) 
3.2 kW 45 KW 
(short duration) , 
100 mm 80 mm 
10—90 mm 0—90 mm 
600 mm 880 mm 
470 mm 380 mm 
(2-410 mm support) 
about 1 m°/h 1 m?/h 
about 230 kg about 225 kg 
(4-23 kg support) 
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Fig. 3b Magnetisation curves of Lange’s electromagnet poles diameter Ø = 80mm . zi d 
weight is the same as that of Lange’s electromagnet and the power required is three 
times less. Clearly, for more specialized purposes, magnets of still better performanc 3 
may be designed. The present one has been conceived for laboratory work, and : 
adapted to the requirements met with in spectroscopy. The supply voltage is standar 
220V. In its present form and throughout. the range of medium field strengtl 
MUR pe is superior. to that of the EE hitherto, in A E 
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